Planet Migration 

To Appear in "Exoplanets," ed. S. Seager, Univ. Arizona Press. 



Stephen H. Lubow 

Space Telescope Science Institute 



Shigeru Ida 

Tokyo Institute of Technology 

Planet migration is the process by which a planet's orbital radius changes in time. The 
main agent for causing gas giant planet migration is the gravitational interaction of the young 
planet with the gaseous disk from which it forms. We describe the migration rates resulting 
from these interactions based on a simple model for disk properties. These migration rates are 
higher than is reasonable for planet survival. We discuss some proposed models for which the 
migration rates are lower. There are major uncertainties in migration rates due to a lack of 
knowledge about the detailed physical properties of disks. We also describe some additional 
forms of migration. 



1. INTRODUCTION 

Planet formation occurs in disks of material that orbit 
young stars (see Chapter 13). At early stages of evolution, 
the disks are largely gaseous and have masses of typically a 
few percent or more of the stellar masses. At later times, af- 
ter several lO^yr, the gas disperses and disks largely consist 
of solid material. Such disks last for timescales of ^ lO^yr. 
The interactions of a young planet with its surrounding disk 
affects the planet's orbital energy and angular momentum. 
One consequence of such interactions is that a planet may 
move radially through the disk (toward or away from the 
central star), a process called migration that is the topic of 
this chapter In this chapter, we mainly concentrate on mi- 
gration caused by gaseous disks. Although gaseous disks 
survive only a small fraction of the lifetime of the star and 
planet (^ 10^°yr), calculations of disk-planet interactions 
described in this chapter show them to be strong enough to 
have caused substantial migration. In fact, a major prob- 
lem is that the effects of migration are typically predicted 
to be too strong to have allowed the formation of gas gi- 
ant planets to complete. That is, the typical timescales for 
migration are found to be shorter than both the lifetimes 
of gaseous disks and the predicted formation timescales of 
gas giant planets in the so-called core accretion model of 
~ lO^yr (see Chapter 13). Migration models applied some 
simple disk structures predict that a forming planet should 
fall into (or very close to) the central star before it grows in 
mass to become a gas giant (Ward, 1997). 

Early theoretical studies of planet migration, such as 
Goldreich & Tremaine, 1980, understandably concentrated 
on the role of migration for the planets in our solar system. 
These studies suggested that there was substantial migration 
of Jupiter while it was immersed in the solar nebula. Migra- 
tion should have occurred both at Jupiter's early stages of 



formation when it was a solid core and at later stages when 
it was a fully developed gas giant planet interacting with 
the solar nebula. The migration timescales were estimated 
to be much shorter than lifetime of the nebula, although the 
direction of migration was not determined. However, the 
evidence for migration was not apparent. The location of 
Jupiter is just outside the snow-line, where conditions for 
rapid core formation are most favorable (see Chapter 13). 
This situation seemed to suggest that Jupiter formed and re- 
mained near its current location. Therefore, a plausible con- 
clusion was that migration did not play an important role for 
Jupiter and therefore perhaps for all planets. 

Evidence for the importance of migration changed dra- 
matically in 1995 with the discovery of the first giant planet, 
51 Peg b, that has an orbital period of only about 4 days 
(Mayor and Queloz, 1995). Subsequent discoveries of 
many giant planets with orbital radii substantially smaller 
than Jupiter's added to the evidence (e.g., Butler et al, 
2006). Their existence suggested that migration had oc- 
curred, since giant planet formation close to a star is not 
likely to occur (Bodenheimer et al, 2000; Ida & Lin, 2004). 

We might crudely think of a planet orbiting in a gaseous 
disk around a star as being similar to a satellite in orbit about 
the Earth. Atmospheric gas drag on the satellite causes it to 
lose angular momentum and spiral down (migrate inwards) 
towards the Earth. However, the disk-planet situation is 
somewhat different. The planet and disk are both in or- 
bit about the star, and there is a relatively small difference 
between the rotational speeds of the planet and its neigh- 
boring gas. Although the velocity difference can produce 
an important level of gas drag for small mass solid objects, 
gas drag does not produce the dominant torque for objects 
as large as planets, as is discussed further in Section 2.1. 
For disk-planet interactions, gravitational torques play the 
dominant role. 
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The gravitational interactions between a planet and a 
gaseous disk are difficult to analyze in a precise way. But 
a general picture has developed. The results depend some- 
what on the disk structure. We consider the disk to be what 
is called an accretion disk. Accretion disks occur in many 
astronomical objects, such as binary star systems and ac- 
tive galactic nuclei. Such disks typically arise as gas flows 
towards a central body with some rotation about that body. 
Due to its angular momentum, the gas cannot fall directly 
onto the central object and instead forms a disk in near cen- 
trifugal balance. Observations of such systems, including 
disks around young stars, reveal evidence that the gas is 
not simply orbiting about the star, but also flowing inward 
(accreting) towards the central object. The inflow velocity 
within the disk is very small compared to the orbital ve- 
locity. But it is sufficient to reveal various observational 
signatures as a consequence of this mass flow towards and 
onto the central object. In particular, the inflow can result in 
the emission of radiation from the accreting gas as it moves 
deeper in the gravitational potential of the central object. 
In the case of interest to planet migration, the gaseous disk 
orbits about a young central star. These so-called protoplan- 
etary disks are described in Chapter 1 1 . 

Turbulence within the disk is thought to be the main 
driver of the accretion. The disk does not rotate about the 
star at a constant rate with distance from the star. Instead, 
as a consequence of centrifugal balance, it undergoes dif- 
ferential rotation, and the gas rotation rate decreases with 
distance from the star. The gas orbits at (nearly) the so- 
called Keplerian orbital frequency at orbital radius r from 
the star of mass Ms given by 

In a simple model, the disk turbulence can be considered 
to cause an effective friction. Neighboring rings of gas 
in the disk interact by frictional forces as they rub against 
each other due to the differential rotation. This friction pro- 
duces orbital energy losses and a torque that causes nearly 
all the material within the disk to move in (accrete), while 
angular momentum is transported outward {Lynden-Bell & 
Pringle, 1974). Instabilities within the disk are the cause 
of the turbulence. The most likely candidates for causing 
the instabilities are disk self-gravity (e.g., Paczynski, 1978; 
Lodato & Rice, 2004) and magnetic fields (Balbus & Haw- 
ley, 1991). Disk self-gravity instabilities are more likely to 
be important at the earliest stages of disk evolution, while it 
is more massive. 

For small mass planets, the density structure of the disk 
is largely unaffected by the presence of the planet. This 
situation leads to the so-called Type I regime of planet mi- 
gration. In this simple model, the disk density structure is 
determined by the disk turbulence. The disk differential ro- 
tation implies that gas that lies outside the orbit of the planet 
moves more slowly than gas inside the orbit of the planet. 
The planet disturbs the disk by its gravity as gas passes near 
it. The disk disturbances act back on the planet and cause 



a net torque to be exerted on the planet. This torque causes 
the planet to migrate. An analysis, such as in Section 2.2, 
shows that the migration timescale (orbit decay timescale) 
for a 10 Earth mass planet embedded in a disk that is some- 
what similar to the minimum mass solar nebula (Hayashi, 
1981) is only about lO^y (see Fig.[l]i. Furthermore, the mi- 
gration timescale varies inversely with planet mass. This 
result has the consequence that planets migrate even faster 
as they grow in mass. The shortness of this timescale poses 
a challenge to our understanding of migration, as discussed 
earlier. 

For a sufficiently high mass planet, its tidal torques can 
overpower the viscous torques associated with disk turbu- 
lence, and the planet can alter the disk density in its vicinity. 
The result is that the planet opens a gap in the disk about its 
orbit. In this situation, we have the so-called Type II regime 
of planet migration. The disk-planet interactions in this case 
are still gravitational. But the gravitational torque on the 
planet is reduced because there is less material close to the 
planet compared to the Type I regime. In the limiting case 
that the planet mass is small compared to the disk mass, the 
planet moves inward like a test particle along with the ac- 
cretion inflow of the gas. The disk, being much more mas- 
sive than the planet, sets the migration rate of the planet to 
its inflow rate. The disk density near the gap region adjusts 
so that the torque on the planet matches the disk's inflow 
speed. 

The mass required for the planet to open a gap, as well 
as the migration rate, depends on the level of disk turbu- 
lence. For typically adopted disk parameters, the transition 
to Type II migration occurs at a few tenths of a Jupiter mass. 
The migration timescales are about lO^yr. This timescale 
is also short compared to the disk lifetimes, but is much 
longer than it would have been if the planet had continued 
to undergo Type I migration (see Fig. [T]). In this idealized 
Type II regime of high disk-to-planet mass ratio, the planet 
migration rate is independent of planet mass, since it is de- 
termined by the disk only. But for less extreme disk-to- 
planet mass ratios, the inertia of the planet plays a role and 
the migration rate decreases somewhat with planet mass. 

The determination of planet migration rates is not easy 
and various aspects remain controversial. There are sev- 
eral reasons for this. Although the interaction between a 
planet and its surrounding disk is gravitational, the inter- 
action behaves in a way that is reminiscent of friction (dy- 
namical friction). Gas that lies somewhat exterior (interior) 
to the planet's orbit rotates slower (faster) than the planet 
and provides an inward (outward) torque on the planet as 
a consequence of the interaction. The net torque therefore 
involves a competition between nearly equal and opposite 
torques. Accurate calculations of the inward and outward 
torques are then required to determine the migration rate 
and the direction of migration. Another complication is that 
the disk-planet interactions generally get stronger closer to 
the planet. The results therefore depend on the details of 
how these torques become limited near the planet. 

Another issue involves the effects of material that lies 
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very close to the planet. This so-called coorbital region in- 
volves gas that does not circulate past the planet (see Fig.|2]). 
Instead, the gas undergoes abrupt changes in its angular mo- 
mentum due to encounters with the planet that take it from 
interior of the planet's orbit radius to exterior of the planet's 
orbital radius and visa-versa. The torques resulting from 
this region are quite different from those involving more 
weakly perturbed material that lies outside this region. The 
two regions, the weakly deflected and strongly deflected re- 
gions, are dynamically quite different and require separate 
descriptions. Yet another issue is that the torques depend 
on the detailed conditions within the disk. The disk struc- 
tural variations, i.e., how the disk density and temperature 
vary in radius, can play an important role in determining the 
torque. The detailed structural properties of disks around 
young stars are not well understood and could be compli- 
cated (e.g., Armitage et al, 2000; Terquem, 2008). The re- 
sulting structures are likely more complicated than the sim- 
ple power-law density distributions of the minimum mass 
solar nebula model. The level of disk turbulence in a pro- 
tostellar disk influences the disk structure and the rate of 
planet migration. The amount turbulence can vary with dis- 
tance from the star and in time, but it is not well understood. 

In Section 2 we describe calculations of migration rates 
due to disk-planet interactions. Section 2.1 discusses the 
effects of gas drag. Section 2.2 describes the migration rate 
of a low mass planet interacting with a disk, modeled sim- 
ply as a set of ballistic particles that undergo mild or weak 
interactions with the planet. This leads to an estimate of 
the migration rate, but the approximations are not accurate 
enough to determine the direction of migration. Sections 
2.3 discusses the torques due to material that comes very 
close to the planet, torques in the coorbital region. Sec- 
tion 2.4 discusses the more accurate determination of Type I 
migration rates for a gaseous disk and the factors affecting 
migration. Section 3 describes some of the outstanding is- 
sues: how planet migration may be slower than estimated in 
Section 2, some effects in planet migration that are missing 
from the idealized models in Section 2, some forms of mi- 
gration not involving a gaseous disk, and some issues about 
the techniques used in carrying out numerical simulations. 
Section 4 describes some future prospects. 

2. Migration Rates Due To Disk-Planet Interactions 

2.1 Aerodynamic Gas Drag 

Aerodynamic gas drag provides the dominant influence 
on the orbital evolution of low mass solids embedded in a 
gaseous disk that orbits a star (Weidenschilling, 1977; Cuzzi 
& Weidenschilling, 2006). The gas is subject to a radial 
pressure force, in addition to the gravitational force of the 
star. For a disk with smooth structural variations in radius, 
this force induces slight departures from Keplerian speeds 
of order {H/rYVlr, where H is the disk thickness and il is 
the angular speed of the disk at radius r from the star. For a 
thin disk {H <C r) whose pressure declines in radius, as is 



typically expected, the pressure force acts radially outward. 
The gas rotation rate required to achieve centrifugal bal- 
ance is then slightly below the Keplerian rate. Sufficiently 
high mass solids are largely dynamically decoupled from 
the gas and orbit at nearly the Keplerian rate, but are sub- 
ject to drag forces from the more slowly rotating gas. The 
drag leads to their orbit decay. Both the drag force and the 
inertia increase with the size of the sohd, R. The drag force 
on the object increases with its area R^), while its iner- 
tia increases with its volume (^ i?'^). In the high mass solid 
regime, the dominance of inertia over drag causes the or- 
bital decay rate to decrease with object size as 1/i? or with 
object mass as 

The orbital migration rate due to gravitational interac- 
tions between an object and a disk increases with the mass 
of the object, as we will see in the next section. There is 
then a cross-over mass above which the orbital changes due 
to disk gravitational forces dominate over those due to drag 
forces. For typical parameters, this value is much less than 
an Earth mass, IM^ {Hourigan & Ward, 1984). Conse- 
quently, for the purposes of planet migration, we will ignore 
the effects of aerodynamic gas drag. 

2.2 Torques Due To Mildly Perturbed Particles 

As an initial description of gravitational disk-planet in- 
teractions, we consider the disk to consist of noninteracting 
particles, each having mass M much less than the mass of 
the planet A/p. The disk is assumed to extend smoothly 
across the orbit of the planet without a gap. The planet is 
then in the Type I regime of planet migration, as discussed 
in the Introduction. The model leads to insights about sev- 
eral issues concerning planet migration. In particular, it ex- 
plains why the migration rate increases with planet mass 
and allows us to estimate the magnitude of the migration 
rate. This estimated migration rate is numerically close to 
what more detailed calculations reveal for simple disk struc- 
tures. The approximations fall short of allowing a determi- 
nation of whether the migration is inward or outward and 
the determination of the detailed dependence of the migra- 
tion rate on disk properties, such as the density and temper- 
ature distributions. The treatment describes material that 
orbits about the star and undergoes a small deflection by 
the planet. It does not describe the effects of material that 
lies very close to the planet in the so-called coorbital re- 
gion, where the deflections are large (see Fig. |2]i. We will 
consider a model for that region in Section 2.3. 

In this model, the particles are considered to encounter 
and pass by the vicinity of a planet that is on a fixed circu- 
lar orbit of radius a about the star of mass Mg. As a result 
of their interactions, the planet and disk exchange energy 
and angular momentum. This situation is a special case of 
the famous three-body problem in celestial mechanics. The 
tidal or Hill radius of the planet where planetary gravita- 
tional forces dominate over stellar and centrifugal forces is 
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where radius i?H is measured from the center of the planet. 

Consider a particle that approaches the planet on a cir- 
cular orbit about the star with orbital radius r sufficiently 
different from a to allow it to freely pass by the planet with 
a small deflection. This condition requires that the clos- 
est approach between the particle and planet ^ |r — a| 
to be somewhat greater than a few times i?H- The solid 
line in Fig. [3] shows the path of a particle deflected by a 
planet whose mass is lO^^Ms (corresponding to 0.3 Earth 
masses for a planet that orbits a solar mass star), Hill ra- 
dius i?H — 0.007a, and orbital separation r — a ~ 3.5i?H- 
The planet and particle orbit counter-clockwise in the iner- 
tial frame with angular speeds fip and ri(r), respectively. 
In the frame of the planet, the particle moves downward in 
the figure, since its angular speed is slower than that of the 
planet {Vl{r) < flp for r > a). 

To estimate the angular momentum change of a particle 
like that in Fig. [3] we consider a Cartesian coordinate sys- 
tem centered on the planet and corotates with it. The x axis 
lies along a line between the star and planet and points away 
from the star. We consider in detail the fate of a particle 
that starts its approach toward the planet on a circular orbit 
of orbital radius r that lies outside the orbit of the planet, 
X — r — a > 0. The planet and the pre-encounter particle 
are on circular Keplerian orbits whose orbital frequencies 
are given by equation ([T]). We later discuss what happens to 
particles that begin their encounter at orbital radii smaller 
than the planet's (x < 0). The dashed line in Fig. [3]traces 
the path that the particle would take in the absence of the 
planet, while the solid line shows the path in the presence 
of the planet. In both cases, the particle paths are generally 
along the negative y direction. The two paths are nearly 
identical prior to the encounter with the planet (for y > 0). 

The particle velocity in the frame that corotates with the 
planet is then approximately given by 

dn 3 

V ~ r(ri(r) — rip)ey ~ xa—By ~ —-flpxey, (3) 

where x = r — a, il{r) is given by equation ([T]), and 
flp = fl{a). Upon interaction with the planet, the par- 
ticle is deflected slightly toward it. The particle of mass 
AI experiences a force Fx = -GMMpx/{x'^ + y{tff^^. 
This force mainly acts over a time t when \y{t)\ < x and 
is Fx ~ -GMMp/x"^. From equation (O, it foUows that 
the encounter time At ^ x/v ^ of order the orbital 

period of the planet, independent of x. 

To proceed, we apply the so-called impulse approxima- 
tion, as employed by Lin & Papaloizou, 1979. The approx- 
imation involves the assumption that the duration of the in- 
teraction is much shorter than the orbital period. Since the 
duration of the encounter is of order the orbital period, the 
approximation is only marginally satisfied. Consequently, 
the expressions we obtain cannot be determined with high 



accuracy in this approximation. They contain the proper de- 
pendences on various physical quantities. But the approxi- 
mation does not lead to the correct dimensionless numerical 
coefficients of proportionality for the angular momentum 
change of a particle. Therefore, we suppress the numeri- 
cal coefficients in the analysis. An exact treatment in the 
limit of weak perturbations exerted by a small mass planet 
is given in Goldreich & Tremaine, 1982. 

As a result of the encounter, the particle with a: > 
acquires an x (radial) velocity 
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The particle is then deflected by an angle 
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(4) 



(5) 



after the encounter (see Fig.|3]l. 

To determine the change in angular momentum of the 
particle, we need to determine its change in velocity along 
the 9 direction, that is the same as the y direction near the 
planet. To determine this velocity change, Avy, we ignore 
the effects of the star during the encounter, and apply con- 
servation of kinetic energy between the start and end of the 
encounter in the frame of the planet. The velocity mag- 
nitude V is then the same before and after the encounter, 
although the direction changes by angle 5. Since the par- 
ticle in Fig. [3] moves in the negative y direction, its pre- 
encounter y velocity is —v and its post-encounter y velocity 
is — u cos 5. We then have that the change of the velocity of 
the particle along the y direction is 



Avy = —V COS 5 



■ V, 



(6) 



where v ~ l.bVlpX is the magnitude of the velocity before 
the encounter (see equation (|3]l). We assume that the per- 
turbation is weak, 5 « 1, and obtain from equations (O, 
©, and © 



(7) 



It then follows that the change in angular momentum of the 
particle is given by 



A J MaAv,, ^ Ma^Vlr 



Mp 
Ms 



(8) 



Fig-Slplots the orbital radius as a function of time for the 
particle plotted in Fig. [3] where t = is the time of clos- 
est approach to the planet. Notice that after the encounter, 
the particle orbit acquires an eccentricity, as seen by the ra- 
dial oscillations, and an increased angular momentum, as 
seen by the mean shift of the radius for the oscillations. Al- 
though the particle is initially deflected inward (negative x 
direction, towards the star), as expected by equation (|4|l, it 
rebounds after the encounter to an increased time-average 
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radius and higher angular momentum, as expected by equa- 
tion (O. 

Fig lSlplots the torque on a particle due to the planet as a 
function of time, where i = is the time of closest approach 
between the particle and the planet. The solid curves are for 
particles that follow their actual paths (similar to the solid 
line in Fig.O. The dashed curves are for particles that are 
made to follow the unperturbed paths (similar to the dashed 
line in Fig. [3]i. The net angular momentum change is the 
time-integrated torque. The torque on a particle along an 
unperturbed path (such as the dashed curve in Fig. |3]l is an- 
tisymmetric in the time t. Consequently, there is no net 
change in angular momentum accumulated along this path. 
The change in angular momentum along the unperturbed 
path from t = — ootoi = Ois linear in the planet mass, 
since the force due to the planet is proportional to its mass. 
The departures from antisymmetry of the torque versus time 
result in the net angular momentum change. These depar- 
tures are a consequence of the path deflection (solid curve 
in Fig.O. 

For the case of the closer encounter plotted in the top 
panel of Fig.|6] the departures from the unperturbed case are 
substantial. But for a slightly larger orbit (bottom panel), 
the torque along the perturbed path differs only slightly 
from the torque along the unperturbed path that integrates 
to zero. This behavior is a consequence of the steep de- 
cline of A J with X in equation dHJ. For x greater than a few 
i?H, the angular momentum change acquired before the en- 
counter is nearly equal and opposite to the angular momen- 
tum change at later times t > 0. The net angular momen- 
tum change is given by equation (|8]l and is quadratic in the 
planet mass. The quadratic dependence on Mp is a conse- 
quence of the deviations in torque between the perturbed 
and the unperturbed paths, since the angular momentum 
change along the unperturbed path is zero. These deviations 
involve the product of the linear dependence of the force on 
planet mass with the linear dependence of the path deflec- 
tion on planet mass (equation (|5]l) and so are quadratic in 
Mp. This quadratic dependence of torque on planet mass is 
important. It leads to the conclusion we obtain later that the 
migration rate of a planet is proportional to its mass. 

The particle in Fig.|3]gains angular momentum, as pre- 
dicted by equation dHJ. The reason is that the path deflec- 
tion occurs mainly after the particle passes the planet, that 
is for y < 0. The deflection takes the particle closer to the 
planet than would be the case along the unperturbed path. 
The planet then pulls the particle toward positive y, caus- 
ing it to gain angular momentum. Just the opposite would 
happen for a particle with r < a, for a; < in equation 
(O. The particle would approach the planet in the positive 
y direction, be deflected towards the planet for y > 0, and 
be pulled by the planet in the negative y direction, caus- 
ing it to lose angular momentum. The interaction behaves 
somewhat like friction. The particle gains (loses) angular 
momentum if it moves slower (faster) than the planet. The 
angular momentum then flows outward as a result of the in- 
teractions. That is, for a decreasing angular velocity with 



radius {dft/dr < 0) as in the Keplerian case, a particle 
whose orbit lies interior to the planet gives angular momen- 
tum to the planet, since the planet has a lower angular speed 
than the particle. The planet in turn gives angular momen- 
tum to a particle whose orbit lies exterior to it. 

The planet in effect pushes material away from its orbit. 
A particle outside the orbit of the planet is forced outward 
(in a time-averaged sense) as it gains angular momentum, 
and a particle inside is forced inward as it loses angular mo- 
mentum. From this point of view, the gravitational effects 
of the planet behave in a repulsive manner. 

Fig. |6] shows the results of numerical tests of equation 
(HJ. It verifies the dependence of A J on x and Mp. Depar- 
tures of the expected dependences (solid lines) occur when 
X ~ 3i?H. At somewhat smaller values of x, the particle 
orbits do not pass smoothly by the planet. Instead, they lie 
within the coorbital region where they periodically undergo 
strong deflections, as seen in Fig. [2] We omit this region 
from current consideration and consider it later in Section 
2.3. 

We now determine the torque on the planet for a set of 
particles that form a continuous disk (of zero thickness) that 
lies in the orbit plane of the planet with surface density 
(mass per unit area) E that we take to be constant in the 
region near the planet. The particle disk provides a flux of 
mass (defined to be positive) past the planet between x and 
X + dx 

dM \vy\i:dx -^^ripxdx. (9) 



We evaluate the torque Tout on the planet due to disk 
material that extends outside the orbit of the planet from 
r — a + Ar to oo or x from Ar to oo, where Ar > 0. 
We use the fact that the torque the planet exerts on the disk 
is equal and opposite to the torque the disk exerts on the 
planet. We then have 



Tout ^ 
Tout = 



AJ dM 

M dx 



dx. 



(10) 
(11) 



where Ct is a dimensionless positive constant of order 
unity and A J is evaluated through equation ([8]l. The torque 
on the planet due to the disk interior to the orbit of the 
planet from x ~ — oo to a; = Ar with Ar < evaluates to 

Tin — —Tout or 



Tu 



= CxSf^pa^ 



Mp 
Ms 



(12) 



The equations of motion for particles subject to only 
gravitational forces are time-reversible. We saw that par- 
ticle in Fig. [3] gains angular momentum from its interac- 
tion with the planet. But, if we time-reverse the particle- 
planet encounter in Fig. |3] we see that the eccentric orbit 
particle would approach the planet (both on clockwise or- 
bits) and then lose angular momentum (apply t — > —t in 
Fig. IHi. What determines whether a disk particle gains or 
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loses angular momentum? We have assumed that the parti- 
cles always approach the perturber on circular orbits. The 
particles periodically encounter the gravitational effects of 
the planet. The closer the particle orbits to the planet, the 
smaller the relative orbital speeds and the longer the time 
between encounters. From equation (|3]l, it follows that the 
time between encounters r ~ 2na/v is estimated as 

where P is the planet orbital period. Since the encounters 
are close |a;| ^ a, the time between encounters is long 
compared to the encounter time ^ P, that is r ^ P. But 
we saw in Fig.|4]that the particles acquire eccentricity after 
the encounter. For this model to be physically consistent, 
we require that this eccentricity damp between encounters 
with the planet. The eccentricity damping produces an ar- 
row of time for the angular momentum exchange process 
that favors circular orbits ahead of the encounter as shown 
in Fig. |3] resulting in a gain of angular momentum for par- 
ticles that lie outside the orbit of the planet. 

Equations ( fTTT i and (fTZt have important consequences. 
The torques on the planet arising from the inner and outer 
disks are quite powerful and oppose each other This does 
not mean that the net torque on the planet is zero because we 
have assumed perfect symmetry across r = a. The symme- 
try is broken by higher order considerations. For this rea- 
son, migration torques are often referred to as differential 
torques. 

Since the torques are singular in Ar, they are domi- 
nated by material that comes close to the planet. Conse- 
quently, the asymmetries occur through differences in phys- 
ical quantities at radial distances A?- from the planet. In a 
gaseous disk, the effects of both temperature and density 
variations in radius can play a role in the asymmetry, as 
well as asymmetries associated with the differential rota- 
tion of the disk. We cannot describe the asymmetries due 
to temperature in this model because temperature is not de- 
scribed by the ballistic particles. For this simple model of 
a (pressure-free) particle disk, we consider for example the 
effect of the density variation in radius that we have ignored 
up to this point, although the other asymmetries are impor- 
tant. It does not matter which quantity is considered for the 
purposes of obtaining a rough expression for the net torque. 

We expand the disk density in a Taylor expansion about 
the orbit of the planet and obtain that 



2Ar— , 

dr 



(14) 
- Ar 



where Ejn and Sout are the surface densities at r = a - 
and r — a + Ar, respectively. Consequently, for \dTj/dr\ ^ 
Sp/ a, we expect that the sum of the inner and outer torques 
to be smaller than their individual values by an amount of 
order Ar/a. Similar considerations apply to variations in 
other quantities. That is, we have that the absolute value of 
the net torque T on the planet is approximately given by 



\Ar 



Ms 



(16) 
(17) 



The above equation for the torque (equation ( fTTj i) must 
break down for small Ar, in order to yield a finite result. 
For values of | Ar| < 3i?H, the particles become trapped in 
closed orbits in the so-called coorbital region (see Fig. |2]l. 
We exclude this region from current considerations, since 
the torque derivation we considered here does not apply in 
this region. In particular, the assumption that particles pass 
by the planet with a small deflection is invalid in this region. 
Using equation (T% with | Ar| ^ i?H, we obtain a torque 



in 



4/3 



(18) 



Another limit on Ar comes about due to gas pressure. 
One effect of gas pressure is to cause the disk to have a 
nonzero thickness H. The disk thickness is measured by its 
scale height out of the orbital plane. The disk scale height 
is determined by the force balance in the direction perpen- 
dicular to the orbit plane, the z or vertical direction. Gas 
pressure forces act to spread the gas in this direction, while 
gravitational forces act to confine it. We consider the ver- 
tical structure of an axisymmetric disk, unperturbed by a 
planet. There are no motions in the vertical direction and 
the disk is said to be in hydrostatic balance. 

The vertical hydrostatic balance condition can be written 

as 

dp{r, z) 



dz 



P9z 



(19) 



where p{r, z) is the gas pressure, p{r, z) is the gas density, 
and gz {r, z) is the vertical gravity. The vertical gravity is the 
z component of the gravitational force per unit disk mass 
due to the star that is equal to 



GMsZ 



-n^ir)z, 



(20) 



where is given by equation ([T]). For a vertically isother- 
mal disk, a disk whose temperature and therefore sound 
speed c(r) is independent of z, we have 



p(r,z) = p(r,z)c^(r). 



(21) 



Substituting equations (|20] | and ( 1211 1 into equation iT% . we 
obtain 



and 



p{r, z) = p(r, 0) exp 
c(r) 



H{r) 



(22) 



(23) 



(15) 



where c is the gas sound speed that is a function of disk 
temperature. Disk thickness is a measure of the impor- 
tance of gas pressure. For conditions in protostellar disks. 
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the disk thickness to radius ratio (disk aspect ratio) is typ- 
ically 0.03 < H/r < 0.1, corresponding to a gas sound 
speed of about 1.5 km s^^ at lAU from a solar mass star. 

For an axisymmetric disk of nonzero thickness, the sur- 
face density S = pdz evaluates to 

/oc 
p(r, z) dz = V2^H{r)p{r, 0), (24) 
-oo 

where p{r, z) is defined by equation (l22T i. The gas den- 
sity is in effect smeared over distance H out of the orbit 
plane. Near the planet, the gas gravitational effects are then 
smoothed over distance H. Distance Ar is in effect lim- 
ited to be no smaller than ^ H. For Ar ^ H, the torque 
expression ([TtI i is then estimated as 

Which form of the torque applies (equation ( fTsT i or (IZST i) 
to a particular system depends on the importance of gas 
pressure. For H < i?H, gas pressure effects are small com- 
pared to the gravitational effects of the planet near the Hill 
radius. Consequently, we expect equation ( fTSl l to be appli- 
cable in the case of relatively weak gas pressure, H < i?H, 
and equation (|25] l to be applicable otherwise. For typically 
expected conditions in gaseous protostellar disks, it turns 
out that equation (IZSl l is the relevant one for planets under- 
going Type I migration. 

The migration rate f-p/a = T/ Jp for a planet with angu- 
lar momentum Jp is then linear in planet mass, since T is 
quadratic while Jp is linear in planet mass. Therefore, the 
Type I migration rate increases with planet mass, as asserted 
in Section 2.1. This somewhat surprising result that more 
massive planets migrate faster is in turn a consequence of 
the quadratic variation of A J with planet mass in equation 
([8]l. This quadratic dependence occurs because the possible 
linear dependence of A J on planet mass vanishes due to the 
antisymmetry of the torque as a function of time along the 
unperturbed particle path, as discussed after equation dHJ. 

Based on equation ( |25] | with typical parameters for the 
minimum mass solar nebula at the location of Jupiter a — 
5AU (E = 150gcm-2, 17 = 1.8 x 10-*s-\ and H = 
0.05a), we estimate the planetary migration timescale Jp/T 
for a planetary core of 10 Earth masses embedded in a mini- 
mum mass solar nebula as 4 x lO^yr. This timescale is short 
compared to the disk lifetime, estimated as several times 
10^ yr, or the Jupiter formation timescale of > 10^ yr in 
the core accretion model. The relative shortness of the mi- 
gration timescale is a major issue for understanding planet 
formation. Since migration is found to be inward for sim- 
ple disk models, as we will see later, the timescale dispar- 
ity suggests that a planetary core will fall into the central 
star before it develops into a gas giant planet. Research 
on planet migration has concentrated on including addi- 
tional effects such as gas pressure and on improving the 
migration rates by means of both analytic theory and multi- 
dimensional simulations. 



A more detailed analysis reveals that the torque in equa- 
tion (|25] | does provide a reasonable estimate for the mag- 
nitude of migration rates in gaseous disks in the so-called 
Type I regime in which a planet does not open a gap in the 
disk. However, the derivation of the torque in this section 
is not precise enough to determine whether the migration is 
inward or outward (i.e., whether T is negative or positive). 
The density asymmetry about r = a (see equation (fT4l i) 
typically involves a higher density at smaller r, as in the 
case of the minimum mass solar nebula. This variation sug- 
gests that torques from the inner disk dominate, implying 
outward migration, as was thought to be the case in early 
studies. But this conclusion is incorrect. As we will dis- 
cuss later, inward migration is typically favored, at least for 
simple disk models. We have not included the effects of 
gas temperature and gas pressure. A disk with gas pres- 
sure propagates density waves launched by the planet. The 
analysis in this section has only considered effects of mate- 
rial that passes by the planet. In addition, there are effects 
from material that lies closer to the orbit of the planet (see 
Fig.|2]l. This region can also provide torques. We have also 
assumed that the disk density is undisturbed by presence of 
the planet. Feedback effects of the disk disturbances and 
gaps in the disk can have an important influence on migra- 
tion. Finally, there are other physical effects such as disk 
turbulence that should be considered. We will consider such 
effects in subsequent sections. 

2.3 Coorbital Torques 

Thus far we have considered torques that arise from gas 
that passes by the planet in the azimuthal direction (the y 
direction in Fig. |3]l with a modest deflection. Gas that re- 
sides closer to the planet, |r — a| < 3i?H, in the so-called 
coorbital region, does not pass by the planet. Instead, it fol- 
lows what are called librating orbits in the corotating frame 
of the planet as seen in Fig. |2] We will consider librating 
orbits of particles that reach close to a planet that lies on a 
fixed circular orbit, the so-called horseshoe orbits. A parti- 
cle at position o+ in Fig.|2]is in circular motion outside the 
orbit of the planet. It moves more slowly than the planet and 
approaches it. The particle in this case is more strongly per- 
turbed by the planet than the more distant particle in Fig. [3] 
It gets pulled inward by the planet, causing it to change to 
a circular orbit interior to the planet's orbit once it reaches 
position i+. Since the particle is now moving faster than 
the planet, it moves away from it. In the process, the par- 
ticle has executed a U-turn. At a later time, the particle 
approaches the planet from behind, at position i— , and the 
planet pulls the particle outward, causing it to make a sec- 
ond U-turn to position o— . The particle is then back to the 
initial outer radius that it started on before both encounters. 
The particle path is closed in the corotating frame of the 
planet. 

With each change in angular momentum of the particle, 
there is an equal and opposite change in angular momentum 
of the planet. But, whatever angular momentum is gained 



7 



by a particle as it changes from the inner radius to outer 
radius Tq is lost when it later encounters the planet and shifts 
from To to ri. The reason is that the particle follows a peri- 
odic orbit in the frame of the planet, so there is no change 
in angular momentum over a complete period of its motion. 
As a result, angular momentum changes do not grow over 
long timescales. At any instant, the largest angular momen- 
tum change possible over any time interval is that acquired 
in the last particle-planet encounter. Therefore, the time- 
averaged torque on the planet due to a particle drops to zero 
over times scales longer than the period of its motion. So al- 
though particles on closed horseshoe orbits can come close 
to the planet, the orbital symmetry limits the torque that 
they exert on the planet. 

For a set of particles, there is a further type of torque 
cancellation that occurs. Recall that the time period be- 
tween particle-planet encounters (the hbration timescale) 
varies inversely with the particle-planet orbital separation 
as 1/x (see equation (fTsli). So even if a set of particles at 
different radii are initially lined up to encounter the planet at 
the same time and produce a large torque, they will eventu- 
ally drift apart in azimuth (phase) and encounter the planet 
at different times. After a while, at a time when one parti- 
cle gains angular momentum another will be losing it. This 
randomization, called "phase mixing", leads to a drop in 
torque over time. Fig.|2]shows the angular momentum evo- 
lution of a set of 60 particles that initially lie on the negative 
X axis within the horseshoe orbit region exterior to the or- 
bit of the planet (r > a). The torque (the derivative of the 
curve) shows both types of cancellations: the time rate of 
change of angular momentum oscillates due to the orbital 
symmetry and the amplitude of the oscillations decreases 
due to phase mixing. 

A space-filling continuous set of small particles will be 
considered to represent coorbital gas. The disk in this case 
is dissipationless and pressureless, and the planet is on a 
fixed circular orbit. As is the case for a set of particles de- 
scribed above, the torque that the gas exerts on the planet 
approaches zero on timescales longer than the time for gas 
to make successive encounters with the planet. 

We now consider the torque caused by gas whose den- 
sity varies with radius with some imposed initial surface 
density distribution S(r). The particle paths in Fig.|2]then 
become streamlines for the gas flow. Gas loses angular mo- 
mentum in going from position o+ to i+. This angular mo- 
mentum is continuously gained by the planet. Similarly, the 
planet continuously loses angular momentum from the gas 
that passes from i_ to o_. The gas spends a relatively short 
time in transition between ri and compared to the time it 
spends between encounters given by equation ( fTSI ), where 
X = a — ri. We consider here the torque exerted on the 
planet on timescales longer than the transition timescale, 
but shorter than the timescale between encounters with the 
planet. 

As gas in the coorbital region approaches from inside 
the orbit of the planet, it carries a mass flux (defined to be 



positive) that we estimate as 



(26) 



where Wi denotes the radial extent (width) of the coor- 
bital horseshoe region interior to the planet's orbital radius. 
Quantities Sj, il;, and ri denote the values of the surface 
density, angular velocity, and orbital radius, respectively, 
at position i~. This position lies at an intermediate radius 
within the horseshoe orbit region and interior to the planet's 
orbit. Gas at this position is executing circular motion, just 
ahead of an encounter with the planet. Similarly, for coor- 
bital gas outside the orbit of the planet, the mass flux is 
estimated by 



riplroWo, 



(27) 



where the subscripts denote the values at a point o+ that lies 
at an intermediate radius within the horseshoe orbit region 
outside the orbit of the planet. 

Due to the orbital geometry, we have that region widths 
satisfy Wi ^ Wq ^ — a ^ a — r^. This gas interacts with 
the planet and undergoes a change in angular momentum 
per unit mass as it flows from position i— to o— 



A J - r^nl 



an^w > 0, 



(28) 



where w — Wo + wi ^ ro — r-i. In going from o+ to i+, 
the gas undergoes an equal and opposite change in angu- 
lar momentum, — AJ. The flows from outside and inside 
the orbit of the planet impart a torque on the planet. This 
torque is equal and opposite to the rate of change of angular 
momentum of gas resulting from the two contributing mass 
fluxes 



Tr. 




(29) 
(30) 



We need to be careful in evaluating the ratio of the mass 
fluxes in the above equation, since departures from unity 
are critical. From equations (|26] | and ( l27b . the ratio of the 
mass fluxes is given by 



(31) 



Ward, 1991 showed that the term in parenthesis on the 
right-hand side is equal to Bo/Bi, where B{r) is the Oort 
constant defined through 2rB = d{r'^Vl)/dr and B{r) = 
0(r)/4 for a Keplerian disk. We then obtain that 



Mi 



Si Bo 
So Bi 



(32) 



Using equation ( l27b . we then approximate Mq appearing in 
the first term on the right-hand side of equation (l30l l as 



Mo ~ Sripaw^. 



(33) 
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(34) 



(35) 



Applying equations (l30l l. (|32] |. and (l33T l. we then obtain an 
expression for the coorbital torque as 

'AE AB\ 

^ a log r 

where AE = So - E; and AB = Bo - Bi. 

Quantity B/E is sometimes called the vortensity, since 
B is half the vorticity, that is half the curl of the unper- 
turbed velocity rileg. The coorbital torque then depends 
on the gradient of the vortensity. Recall that in this situa- 
tion, the torque is not permanent. It oscillates and decays 
to zero, due to phase mixing, as discussed above. The coor- 
bital torque drops to zero or is said to be saturated on long 
timescales. The gas density in the coorbital region becomes 
modified over time so that the vortensity becomes constant 
and the torque approaches zero. 

Gas that flows on the closed streamlines shown in Fig. |2] 
can impart a change in angular momentum on the planet 
that is at most about equal to the gas angular momen- 
tum of the region times its fractional width, w/a (that is, 
Eripfl^tu^). This angular momentum change is due to the 
rearrangement of gas within the horseshoe orbit region from 
some initial state. In this model, gas is trapped within the 
coorbital region and cannot exchange angular momentum 
with the large reservoir in the remainder of the disk. Main- 
taining a coorbital torque over long timescales requires the 
action of a process that breaks the symmetry of the stream- 
lines and permits exchange with the remainder of the disk. 
One such process involves the disk turbulent viscosity. The 
disk viscosity acts to establish a characteristic density distri- 
bution, as discussed in the Introduction. Sufficiently strong 
turbulence can overcome the effects of phase mixing that 
saturate the torque. The torque is then given by equa- 
tion (|35T l. where the density distribution E(7') is determined 
by turbulent viscosity in the viscous accretion disk. The 
vortensity gradient is then maintained at a nonzero value. 
The angular momentum changes within the coorbital re- 
gion are transfered to the remainder of the disk through 
torques associated with frictional stresses caused by disk 
turbulent viscosity (viscous torques). In this way, a steady 
state torque can be exerted on the planet over timescales 
much longer than the libration timescale ^ aP/w. 

We compare the coorbital torque with the net torque due 
to particles outside the coorbital region in a pressureless 
disk, as we discussed in Section 2.2. Taking w ^ i?H in 
equation ( [35l l, we see that the unsaturated coorbital torque 
(taking |c?log (i3/E)/c?logr| ~ 1 and nonzero) is com- 
parable to the torque due to particles outside the coorbital 
region in equation ( fTsT l. It can be shown that this is also 
true for a disk where pressure effects are important, where 
H > i?H- For such a disk, the coorbital torque is generally 
of the same order as the net torque from the region outside 
the coorbital zone, equation ( |25] ). The direction of the coor- 
bital torque contribution depends on the sign of the vorten- 
sity gradient. For the minimum mass solar nebula model. 



E c>c r^^/'^ and B oc r^^^^. Consequently, the coorbital 
torque is zero in that case. With more general density dis- 
tributions of the form 



(36) 



and for smaller values of (3 < 3/2, the coorbital torque on 
the planet is positive. 

We have discussed how the disk turbulent viscosity can 
cause a torque to be exerted in the coorbital region over 
long timescales. But if the amount of viscosity is too small, 
the phase mixing process discussed above can dominate and 
cause the torque to saturate. For the turbulent viscosity to 
prevent torque saturation, the timescale for turbulence to af- 
fect the density in the coorbital region needs to be shorter 
than the hbration timescale on which the torque would drop 
in a nonviscous disk. This criterion imposes a constraint on 
the level of viscosity required. So although they are poten- 
tially powerful, coorbital torques are not guaranteed to play 
a role in planet migration in all cases. 

Although the coorbital torque is of the same order as the 
net torque due to material outside the orbital zone, more 
detailed calculations as described in the next section show 
it typically does not dominate the total torque, at least for 
simple disk conditions. 

Although the description here is self-consistent, coor- 
bital torques are less well understood than the torques in- 
volving noncoorbital material we considered in Section 2.2. 
They are more complicated because they are subject to satu- 
ration. The asymmetry in the coorbital region that prevents 
torque saturation can in principle be produced by effects 
other than turbulence, such as the migration of the planet it- 
self (see Section 3.2). Furthermore, recent simulations sug- 
gest that the thermodynamic state of the gas, not considered 
in this model, can modify the coorbital torque and possibly 
result in outward migration, as is discussed in Section 3.2. 

2.4 Type I Migration 

The torques arising from gas that lies somewhat inside 
the orbit of the planet r < a (but outside the coorbital re- 
gion) act to cause outward migration, while torques arising 
from gas outside the orbit of the planet act to cause inward 
migration, as discussed in Section 2.2. The net migration 
torque was estimated in equation ( l25T l. But due to the in- 
accuracy of the approximations, the sign of the torque and 
its dependence on gas properties, such as the density and 
temperature distributions, could not be determined. More 
accurate analytic calculations of the torque involve solving 
the fluid equations for gas subject to the gravitational per- 
turbing effects of the planet. Both the weakly perturbed 
region (Section 2.2) and coorbital region (Section 2.3) are 
analyzed. The calculations described in this section assume 
a simple disk model. By a simple model we mean that the 
disk has smooth disk density distribution, in the absence of 
the planet, and a smooth and fixed temperature distribution, 
even in the presence of the planet. 
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Consider a cylindrical coordinate system centered on the 
star (r, 9, z). In these calculations, the perturbing potential 
of a circular orbit planet is expanded in a Fourier series as 

$(r, e,t)^Y. *'«(^' «) COS - ^pt)l (37) 

where m is a nonnegative integer, and ilp is the orbital fre- 
quency of the planet. Such an expansion is possible because 
the potential is periodic in azimuth and periodic in time with 
the planet's orbital period. In addition, at fixed radius r (and 
fixed a), the potential depends only on the relative azimuth 
between that of a point and that of the planet Vl-pt, that is 
9 — ilpt. For each azimuthal wavenumber m, the gas re- 
sponse is calculated by means of linear theory. There is a 
torque associated with each m-value due to the effects of the 
density perturbation for each m. The sum of these torques 
determine the net torque on the planet. 

The results show that the gas response is dominated by 
the effects of resonances {Goldreich & Tremaine, 1979, 
1980). There are two types of resonances that emerge: 
Lindblad resonances and corotational resonances. They oc- 
cur in the regions described in Sections 2.2 and 2.3, respec- 
tively. Such resonances also arise in the theory of the spiral 
structure of galaxies and planetary rings. 

The gas response to a particular Fourier potential com- 
ponent m in equation ( l37b is strong at particular radii r„i 
where the Lindblad resonance condition is satisfied. The 
Lindblad resonances occur for gas that periodically passes 
by the planet. They correspond to the so-called mean mo- 
tion resonances of particles in celestial mechanics. For a 
given m-value, the forcing frequency experienced by a par- 
ticle is the time derivative of the argument of the cos func- 
tion in equation ( |37] | along a particle's path 9{t). The forc- 
ing frequency is then given by m{n{r) ~ flp). Consider a 
particle in a circular orbit about a central mass. If the par- 
ticle is momentarily slightly perturbed, it oscillates radially 
at a frequency called the epicyclic frequency, often denoted 
by K{r). This frequency is the free oscillation frequency of 
the particle, analogous to the free oscillation frequency of 
a spring in a simple harmonic oscillator For a Keplerian 
disk, the radial frequency K(r) is equal to the circular fre- 
quency fl{r). (This is why noncircular Keplerian orbits are 
closed ellipses in the inertial frame.) Lindblad resonances 
occur wherever the absolute value of the forcing frequency 
matches the free oscillation frequency. For a Keplerian disk, 
the Lindblad resonance condition, \m{^l — flp)\ — k, sim- 
plifies to 

n{r„.) = (38) 

where ri(r) is given by equation 

For each m-value there are two Lindblad resonances. An 
outer Lindblad resonance (OLR) occurs outside the orbital 
radius of the planet, for which the plus sign is taken in the 
denominator of equation ( l38T l and an inner Lindblad reso- 
nance (ILR), for which the minus sign is taken (see Fig.|8]l. 
At each Lindblad resonance, spiral waves are launched that 



propagate away from the resonance (see Fig. |9]l. These 
waves are similar to acoustic waves, but are modified by 
the disk rotation. One can associate an angular momentum 
and energy with the waves. They carry energy and angu- 
lar momentum from the planet, resulting in a torque on the 
planet. Various processes cause these waves to damp as 
they propagate. As the waves damp, their angular momen- 
tum is transfered to the disk. This situation is reminiscent 
of the case of ocean waves. They are generated by wind 
far from land, but undergo final decay when they break at 
the shore and can cause irreversible changes. Similarly, the 
waves in disks can modify the underlying disk density dis- 
tribution I](r) and open gaps as they damp, although the 
disk turbulence can wash out these effects. 

At the corotation resonance in a Keplerian disk, the con- 
dition is simply that 

^{r-m) = f^p, (39) 

or = a for all m. That is, the corotation resonance oc- 
curs at the orbital radius of the planet. It lies in the region 
of orbits shown in Fig. IH The corotation resonance also 
creates disturbances in the disk, but the disturbances do not 
propagate. Instead they are evanescent and remain trapped 
within a radial region whose size is of order the disk thick- 
ness H. 

The analytic calculations determine a torque for each m- 
value at the inner and outer Lindblad resonances and the 
corotation resonance. The strength of the Lindblad torques 
increases for resonances that lie closer to the planet. The 
closer resonances occur for higher m-values. But due to ef- 
fects of pressure, the resonance condition (l38l l breaks down 
close to the planet and the torque reaches a maximum value 
for mcr ^ r/H where the resonance is at a distance of or- 
der the disk thickness from the planet (see Fig. [8]). Three di- 
mensional effects also weaken the torque close to the planet. 
The reduction arises because the gas is spread over the disk 
thickness H (as was discussed in Section 3.2). The critical 
m-value of met is sometimes called the torque cutoff. This 
reduction of torque for m > mcr impUes that the sum of the 
torques over all m-values remains finite. In this way, these 
calculations avoid the the singularity that we encountered 
in equation (fTTI i involving Ar. 

The results of these calculations show that planets have 
a definite tendency to migrate inward (Ward, 1997). For a 
given m-value, the torque is stronger at the OLR than the 
ILR in a disk, even if we ignore the asymmetries caused 
by density and temperature variations. There are several 
contributions to this density- and temperature-independent 
effect. One such contribution can be seen by noticing that 
for a given m, the OLRs lie slightly closer to the planet's 
orbital radius than ILRs in equation (l38T l. This asymme- 
try, as well as other density- and temperature-independent 
effects, enhances the inward migration effects caused by 
OLRs over the outward migration effects that are caused 
by ILRs. With simple disk structures (smooth disk density 
distributions and smooth, fixed temperature distributions). 
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such effects often dominate the overall torque on the planet 
for typical density and temperature distributions. 

Density and temperature effects can also play an impor- 
tant role in determining the Lindblad torque. The surface 
density distribution in radius influences the torque balance 
in at least two ways. It affects the amount of gas at each 
resonance and it affects the radial pressure force. These 
two effects typically act in opposition to each other and par- 
tially cancel. A declining surface density distribution with 
radius, as is typical, places more gas at the ILRs than the 
OLRs. This effect provides an outward torque contribution, 
since the ILR torques are enhanced. Pressure effects typ- 
ically cause the gas at orbital radius close to a to rotate a 
little more slowly than the planet, as discussed in Section 
2. 1 . Although the primary effect here is gravitational and 
not gas drag as in Section 2.1, the gravitational torque on 
the planet again produces a qualitatively similar effect as 
drag and acts to cause inward migration. This effect can 
also be understood in terms of a slightly decreased angular 
velocity 51 (r) (below the Keplerian rate) in equation ( l38b 
that causes both the ILRs and OLRs to move inward. The 
ILRs are then shifted away from the planet, while the OLRs 
are shifted towards it. The effect of pressure is to favor the 
OLRs, which contribute to inward migration. 

The temperature distribution in radius also has multiple 
effects on the Lindblad torque balance. It affects the gas 
radial pressure force in a similar way that the density dis- 
tribution does, as just described. The pressure effects of 
a radially declining temperature then contribute to inward 
migration. Temperature variations also modify the torque 
cutoff nici ^ a/H, since disk thickness H depends on the 
gas sound speed and therefore temperature (see equation 
( |23]) ). The OLRs, being at a lower temperature (smaller 
H) than the ILRs, have a higher cuttoff rrici. This effect 
favors the OLRs and its inward torques. Both effects of a 
declining temperature with radius, also generally expected, 
contributes to inward migration. 

Detailed 3D linear analytic calculations of the Type 1 mi- 
gration rates have been carried out by Tanaka et al, 2002. 
They assumed that the gas sound speed is strictly constant 
in radius. That is, the gas temperature is assumed to be un- 
affected by the planetary perturbations. For the case of sat- 
urated (zero) coorbital corotation torques, where only dif- 
ferential Lindblad torques are involved, the torque on the 
planet is given by 

T = - (2.34 - 0.10 13) S(a) n{af a' ' (-|) ' , 

(40) 

where /3 is given by equation (|36]|FI The torque on the planet 
resulting from the action of both Lindblad and (unsaturated) 



""^A higher sensitivity to density gradients (more negative coefficient of /3) 
was found in the analysis by Menou & Goodman, 2004. 



coorbital corotation resonances is given by 

T = - (1.36 + 0.54 /3) S(a) niaf (^^) ' (-|) ' . 

(41) 

These migration rates are consistent with the estimate in 
equation ( |25] l. Numerical values for migration timescales 
based on equation dTTT l are plotted in Fig.[T] As mentioned 
above, temperature variations can also contribute to planet 
migration, but are not included in this treatment. Other re- 
sults (such as the 2D calculation by Ward, 1997) suggest 
that the temperature variation should provide an additional 
term within the first parenthesis on the right-hand sides of 
equations ( |40| | and i4l[ that is ^ x, where 

T{r) cx r-^ (42) 

is the temperature distribution with radius (not to be con- 
fused with torque). We typically expect in a disk that 
< ^ < 1.5 and that 0.5 < x < 0.75. Consequently, 
the migration is predicted to be inward and of order the rate 
that would be obtained even if density and temperature gra- 
dients {(3 and x) were ignored. But, such gradients may take 
on different values than assumed here and potentially play 
an important role in slowing migration, as will be discussed 
in Section 3.1. In addition, some further effects, such as 
gas entropy gradients, can modify the coorbital torque. Re- 
cent studies have suggested that such thermal effects could 
change Type 1 migration and may even lead to outward mi- 
gration (Section 3.1; Paardekooper & Mellema, 2006). 

In comparing the saturated coorbital torque (equation 
(|40] |) to the unsaturated coorbital torque (equation i4l[ ). we 
see that the effects of the coorbital torque reduce the inward 
migration rate for /3 < 1.5 and nearly vanish for (3 — 1.5. 
This behavior is consistent with equation ( l35T l, using the 
fact that B cx r^^^. The coorbital torque is then positive 
for (3 < 1.5 and is zero for /3 — 1.5. 

Nonlinear 3D hydrodynamical calculations have been 
carried out to test the migration rates, under similar disk 
conditions (in particular, local isothermality) used to derive 
the analytic model. Figs.fTlandfTOlshow that the migration 
rates agree well with the expectations of the theory. 

We examine the comparison between simulations and 
theory in more detail by comparing torque distributions in 
the disk as a function of disk radius. We define the distribu- 
tion of torque on the planet per unit disk mass as a function 
of radius asdT/dM{r) = 1/ {2TTrY.{r)) dT / dr{r). Fig.E] 
plots dT/dM (scaled as indicated in the figure) as a func- 
tion of radial distance from the planet based on 3D simu- 
lations. The distributions show that the region interior (ex- 
terior) to the planet provides a positive (negative) torque 
on the planet, as predicted in equations (fTTl i and (fT2l i for 
a particle disk. Also, the integrated total torque is nega- 
tive, implying inward migration. The theory predicts that 
the torque density peak and trough occur at distance from 
the planet r — a = Ar ^ tH, where the torque cutoff 
takes effect. For the case plotted in Fig. [TT] that adopts 
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H = 0.05r, the predicted locations agree well with the lo- 
cations of the peaks and troughs in the figure. Furthermore, 
the torque and torque density should scale with the square 
of the planet mass. Although the vortensity gradient is not 
small, since S /B ^ r, we do not find large contributions to 
the torque density from the coorbital region, as is also ex- 
pected. As seen in Fig.[TTl this expectation is well met for 
the two cases plotted. 

2.5 Type II migration 

The planet's tidal torques cause material interior to the 
orbit of the planet (not in the coorbital region) to lose an- 
gular momentum and material exterior to the orbit of the 
planet to gain angular momentum, as we saw in Section 
2.2. The torques then act to clear a gap about the orbit of 
the planet. As discussed in the Introduction, Type II migra- 
tion occurs when the planet mass is sufficiently large that 
tidal forces cause a gap to clear about the orbit of the planet 
{Lin & Papaloizou, 1986). The tidal torques on the disk 
interior or exterior to the planet are estimated by equation 
( fTTT l with Ar ^ H. The tidal torque on the disk that acts to 
open the gap is then estimated as 



To 



Ms 



(43) 



Turbulent viscosity acts to close the gap. The effects of 
turbulence are often described in terms of a kinematic vis- 
cosity I' that is parameterized in the so-called a-disk model 
as 

ly = acH, (44) 

where a is a dimensionless number, c is the gas sound 
speed, and H is the disk thickness, described by equation 
( l23T l {Shakura & Sunyaev, 1973). The value of the key quan- 
tity a is expected to be less than unity because we presume 
that the turbulent motions are slower than sonic and that the 
characteristic length scales for the turbulence are smaller 
than the disk thickness. 

The turbulent viscosity provides a torque to close the gap 
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where Mg is the mass supplied by the disk in closing the 
gap, AJg is the change in the disk angular momentum per 
unit mass required to close the gap, and tg is the time for 
gap closing. We estimate these quantities as 
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where w is the radial extent of the gap, S is the disk den- 
sity just outside the gap, and v is the disk kinematic tur- 
bulent viscosity. The gap closing timescale estimate tg is 



based on the idea that the viscosity acts on a radial diffu- 
sion timescale across a distance w, which then implies the 
quadratic dependence on w. Applying this last set of rela- 
tions to equation (l45T l, we obtain that 



(49) 



independent of w. 

In order to open a gap in the disk, the gap opening torque 
must be greater than the gap closing torque. To > Tc. We 
then obtain a condition on the planet mass required to open 
a gap 
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where Cg is a dimensionless number of order unity. Lin & 
Papaloizou, 1986 estimated that Cg — 2\/T0. For disk pa- 
rameters a — 0.004 and H/r = 0.05, the predicted gap 
opening at the orbit of Jupiter about a solar mass star oc- 
curs for planets having a mass Mp > 0.2Mj, where Afj is 
Jupiter's mass. This prediction is in good agreement with 
the results of 3D numerical simulations (see Fig.fT2l). 

In addition to the above viscous condition, an auxiliary 
condition for gap opening has been suggested based on the 
stability of a gap. This condition is to preclude gaps for 
which steep density gradients would cause an instability 
that prevents gap opening. This condition, called the ther- 
mal condition, is given by the requirement that the Hill ra- 
dius, given by equation (|2]i, be larger than the disk thick- 
ness, i?H ^ H (Lin & Papaloizou, 1986). The critical mass 
for gap opening by this condition is given by 



M^'^ \a 



(51) 



For H/r = 0.05, this condition requires a larger planet 
mass for gap opening than equation ( fSOl l for a < 0.01. A 
condition that combines both equations dSOb and dSTT i has 
been proposed by Crida et al, 2006. Even if both the ther- 
mal and viscous conditions are satisfied and i?H ^ H, a 
substantial gas flow may occur though the gap and onto the 
planet (Artymowicz & Lubow, 1996; Lubow & D'Angelo, 
2006). 

The migration rate of a planet embedded in a gap is quite 
different from the Type I case that we have already consid- 
ered. A planet that opens a gap in a massive disk, a disk 
whose mass is much greater than the planet's mass, would 
be expected to move inward, pushed along with the disk 
accretion inflow. The planet simply communicates the vis- 
cous torques across the gap by means of tidal torques that 
balance them. The Type II migration timescale is then of 
order the disk viscous timescale 



a 

V 



acH 



(52) 



which is ~ 10^ years at Jupiter's orbital radius about the 
Sun for a = 0.004, H = 0.05a, and Vlp = 2tt/\2 yr-^. 
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Therefore, the migration timescale can be much longer 
than the Type 1 migration timescale for higher mass plan- 
ets that open gaps, as is found in simulations (Alp > 
O.lA/j in Figs. [Tl and [TOb. However, this timescale is still 
shorter than the observationally inferred global disk deple- 
tion timescales ^ 10^ — 10^ years. The actual migration 
rate may need to be somewhat smaller, in order to explain 
the abundant population of observed extrasolar giant plan- 
ets beyond lAU (Ida & Lin 2008). 

In practice, the conditions for pure Type 11 migration 
are unlikely to be satisfied. The disk mass may not be 
very large compared to the planet mass and the disk gap 
may not be fully clear of material. However, simulations 
have shown that to within factors of a few, the migration 
timescale matches the viscous timescale, ivis- This result 
holds over a wide range of parameters, provided that the 
tidal clearing is substantial and the disk mass is at least com- 
parable to the planet mass (e.g., Fig.fTOli. 

Another way of understanding the Type II torques is to 
recognize that the distribution of tidal torques per unit disk 
mass, as seen in Fig.[TT]still approximately applies, even if 
the disk has a deep gap that is not completely clear of mate- 
rial. The disk density through the gap region adjusts so that 
the planet migration rate is compatible with the evolution of 
the disk-planet system. 

3. OUTSTANDING QUESTIONS 

We describe some of the major issues involved with the 
theory of planet migration. The topics, descriptions, and 
references are by no means intended to be complete. The 
purpose is to introduce a few of the questions and some of 
the suggested solutions. 

Planet migration is difficult to calculate for a variety 
of reasons, as described near the end of the Introduction. 
There are technical challenges in determining the prop- 
erties of migration torques. For example, the coorbital 
torques are subject to saturation by nonlinear feedback ef- 
fects. Lindblad torques grow in strength with proximity 
to the planet. Their contributions depend on the details of 
how the torques are limited near the planet. For the simple 
disk structures (e.g., moderate turbulent viscosity, fixed and 
smooth temperature distributions, smooth density distribu- 
tions - apart from the density perturbations and gaps pro- 
duced by the planet), the theory has been verified by vari- 
ous multi-dimensional nonlinear hydrodynamic simulations 
(e.g., Figs.[T][TOl andfTTTl. But such simple disk structures 
may not arise in real systems. Indeed, the problems with 
shortness of planet migration timescales suggests that more 
complicated disk structures and physical processes (such as 
more complicated thermal effects) may be important. The 
disk structures influence migration, since migration torques 
are sensitive to gradients of disk properties. But we do not 
have a complete theory for these disk structures. In addi- 
tion, migration depends somewhat on the detailed proper- 
ties of disk turbulence that are not well understood. 

3.1 Limiting/Reversing Type I Migration 



As seen in Fig. [T] the timescales for Type I migration, 
based on the model described in Section 2, are short com- 
pared to gas disk lifetimes of several million years. They 
become even shorter for a disk having a mass greater than 
the minimum mass solar nebula. As emphasized earlier, the 
shortness of the timescales for Type I migration is a serious 
problem for understanding how planets form and survive 
in a gaseous disk. To be consistent with the ubiquity of 
extrasolar gas giants and formation of Jupiter and Saturn, 
some studies suggest that the inward Type I migration rates 
must be reduced by more than a factor of 10 (Alibert et al. 
2005; Rice & Annitage 2005; Ida & Lin 2008). One major 
question is whether there are processes that could slow the 
migration. Several ideas have been proposed. Major un- 
certainties about them arise due to our lack of knowledge 
about the detailed structure of the disks. 

We briefly discuss below a few of the several suggested 
mechanisms for slowing Type I migration. We also discuss 
some processes that may reverse migration, resulting in out- 
ward migration. 

The fundamental tendency for Type I to be rapid can be 
difficult to cancel without artificially fine tuning some other 
effect that provides outward torques. In general, processes 
that provide outward torques tend to result in either rapid in- 
ward or rapid outward migration. However, there are effects 
that provide an opposing torque that increases as the planet 
continues migrating (via feedback effects or "traps"). Such 
processes can naturally maintain halted migration once the 
inward and outward torques are in balance. 

Weak turbulence 

For an inwardly migrating planet, there is a dynami- 
cal feedback effect associated with the radial motion of 
the planet that raises the axisymmetric gas surface density 
within a region interior to the orbit of the planet, and low- 
ers the density within a region exterior to the orbit (Ward, 
1997), in the same sense as a plow operating in the radial 
direction. The planet pushes material away from its orbital 
radius as it migrates. This leads to a pile-up of material 
in the direction of its radial motion (see Fig. [TsT l. The 
feedback then enhances the positive torques that arise in 
the inner disk and slows inward migration. For the feed- 
back to exist, the effects of disk turbulence must be suffi- 
ciently weak, otherwise the density perturbations are erased 
by the effects of turbulent diffusion. This feedback grows 
with planet mass. Above some critical planet mass A/cr, 
the feedback becomes strong enough that the planet can no 
longer migrate. When migration is halted, the density per- 
turbations are initially mild and there is no substantial gap 
(unlike the Type II migration case). There is just a sufficient 
density asymmetry across the orbit of the planet to change 
the competition between inward and outward torques away 
from favoring the inward torques. In subsequent evolution, 
the planet begins forming a gap. 

The value of the critical planet mass depends on 
several factors such as how far from the planet the den- 
sity perturbations occur For disks with H/r ~ 0.05 and 
low turbulent viscosity a < 10^"*, the density perturbations 
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produced by shocked density waves in 2D disks cause mi- 
gration to be sharply reduced at values of Afci ^ lOAf^ (10 
Earth masses) (Rafikov, 2002; Li et al, 2009; see Fig. [T4l) . 
For disk turbulent viscosity parameter a > 10^^, Type I 
migration proceeds with little reduction, since the density 
perturbations are washed out by the effects of disk turbu- 
lence. We have little direct information on what the a 
value should be. Some indirect evidence based on observed 
accretion rates onto young stars suggests that the average 
a > 10^'*. Theory suggests that there may be disk regions 
of low a in magnetically unstable disks (the so-called dead 
zones of low ionization Gammie, 1996). 

Disk Density/Temperature Variations 

The torque on a planet depends somewhat on differences 
in properties of the disk across the planet's orbit, as we have 
seen in Section 2. If the density and temperature smoothly 
decreases in radius at rates typically expected, the outcome 
is inward migration (see equation (HTt and Fig. [U. We ex- 
pect this to be generally true across the disk, but there could 
be regions where the density abruptly changes or where 
the radial temperature gradient is less negative or is even 
positive. To slow migration, the effects of these gradients 
would need to counteract the general trend of inward mi- 
gration discussed earlier that is due to other asymmetric ef- 
fects, such as the asymmetries associated with Keplerian 
rotation. It is possible that rapid changes in disk properties 
with radius could occur as a consequence of strong radial 
variations in disk opacity or turbulent viscosity. They can in 
turn modify the typically adopted gradients of density and 
temperature that led to our estimated migration timescales. 
A planet could then experience slowed migration or even 
be trapped with no further inward migration (e.g., Menou 
& Goodman, 2004; Masset et al, 2006, Matsumura et al, 
2007). However, our current knowledge about disk struc- 
ture is limited. 

Disk temperature variations can also be caused by the 
planet itself. The planet reduces the disk thickness H of 
gas close to it, as a consequence of its gravity. As a result, 
gas close to and just outside the orbit of the planet is further 
exposed to the stellar radiation and experiences additional 
heating by the star Gas just inside the orbit of the planet 
experiences some shadowing. This effect acts to decrease 
the temperature gradient x defined in equation (|42] |. Calcu- 
lations by Jang-Condell & Sasselov, 2005 have shown that 
this effect can reduce the Type 1 migration rates by up to 
about a factor of 2. 

Turbulent Fluctuations 

The effects of disk turbulence are typically described by 
means of a turbulent viscosity, e.g., equation ( l44l i. The 
numerical simulations used in Fig. [T] applied this model. 
Viscosity provides an approximate description of the dy- 
namical effects of turbulent fluctuations that are averaged 
over the small space and short time scales characteristic of 
the turbulence. But, the time-dependent, small-scale den- 
sity fluctuations can give rise to fluctuating random torques 
on the planet that are not described by viscosity. Unlike 
the Type 1 torque that acts continuously in the same di- 



rection, the fluctuating torque undergoes changes in direc- 
tion on timescales characteristic of the turbulence that are 
short compared with the migration timescale. The fluctu- 
ating torque causes the planet to undergo something like a 
random walk. For the effects of the random walk to com- 
pletely dominate over Type 1 migration, the amplitude of 
the fluctuating torque must be much larger than the Type 1 
torque that acts steadily. 

The change in the angular momentum of the planet 
due to a random torque Tr over some time t is given by 
^ -\/]VTRiR , where N = t/tR is the number of fluctuations 
felt by the planet and t-R, is the characteristic timescale for 
the torque fluctuation, perhaps of order the orbital period. 
The planet experiences many torque fluctuations over the 
migration timescale, ^ 1. The change in the planet an- 
gular momentum by the Type 1 torque T is given by ^ Tt. 
For random torques to dominate, we require Tr > \/NT. 
We take t to be the Type 1 timescale of ~ 10^ yr and 
tR ~ 10 yr, so that A^ ^ 10* and we require Tr to be 
> IGGT. Torque Tr depends linearly on the planet mass, 
while the Type 1 torque T increases quadratically with the 
planet mass and the Type 1 migration time t decreases with 
the inverse of the planet mass. Therefore, the random torque 
is then more important for lower mass planets. The nature 
of the random torque depends on the properties of the disk 
turbulence, in particular the amplitude Tr and correlation 
timescale tR for the fluctuating torques, which are currently 
not well determined. 

Some global simulations and analytic models suggest 
that turbulent fluctuations arising from a magnetic insta- 
bility (the magneto-rotational instability Balbus & Hawley, 
1991) are important for migration of lower mass planets 
{Nelson, 2005; Johnson et al, 2006; see Fig. [TSl l. Some 
other simulations, those done for a small region of the 
disk called "box" simulations, suggest that these fluctuat- 
ing torques are not effective enough to play an important 
role in planet migration {Yang et al., 2009). Both the global 
and box simulations are very computationally demanding 
and the current results cannot be regarded as definitive. The 
box simulations have higher resolution than the global ones, 
since they cover a smaller region. However, they may miss 
important effects, if they occur on larger scales than the box 
size. The direct global simulations, such as those in Fig.fTSl 
cover much less time than the Type 1 migration timescale 
for a low mass mass planet and consequently may not re- 
veal the competing effects of unidirectional Type 1 torques. 

Unlike Type 1 migration, the survival of a planet within a 
disk where fluctuating torques dominate is described statis- 
tically. Over time, there is a smaller probability that planets 
survive in the disk. But there is always a nonzero probabil- 
ity of survival. 

However, if turbulent fluctuations are important for mi- 
gration, then the eccentricities of planetesimals are pumped 
up so highly that collisions between them may result in their 
destruction rather than accretion {Ida et al, 2008). There- 
fore, although the turbulent fluctuations may inhibit the in- 
fall of planetary cores into the central star by migration. 
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they may inhibit the build up of the cores necessary for gi- 
ant planet formation in the core accretion model. 
Ordered Magnetic Field 

An ordered magnetic field can affect the nature of disk- 
planet interactions. Terquem, 2003 analyzed the effects of 
an ordered (nonrandom) torroidal magnetic field on Type 1 
planet migration. The magnetic field introduces additional 
resonances, called magnetic resonances, that lie closer to 
the planet than the Lindblad resonances, whose locations 
are given by equation dSSl l. As a result, the magnetic res- 
onances can be stronger than the Lindblad ones. As in the 
case of the Lindblad resonances, a magnetic resonance lo- 
cated interior (exterior) to the orbit of the planet causes 
outward (inward) planet migration. If the magnetic field 
strength falls off sufficiently fast with distance from the 
star, the inner magnetic resonance dominates, and slowed 
or even outward migration can occur For example, in the 
case that the magnetic energy density is comparable to the 
gas thermal energy density in the disk, an outward torque on 
the planet can be produced, if the magnetic-to-gas energy 
density ratio falls off faster than r^^. How this ratio this 
varies within a disk is not known, although such a variation 
is quite plausible in some situations. If the magnetic field is 
responsible for the disk turbulence, then the torroidal field 
would contain a fluctuating component that complicates the 
outcome. 

Migration Driven By Nonisothermal Effects in the Coor- 
bital Region 

The disk is heated by the central star and by viscous 
turbulent dissipation. Many studies of disk-planet interac- 
tions simplify the disk temperature structure to be locally 
isothermal. In this approximation the temperature distribu- 
tion depends on radius in a fixed, prescribed manner. Such 
a situation can arise if the optical depth of the disk is low 
enough that it efficiently radiates any excess energy due to 
compression caused by interactions with the planet. The 
locally isothermal assumption is frequently applied in nu- 
merical simulations and was applied in obtaining equation 
(HTI) . The behavior in the isothermal limit tends to suggest 
that coorbital torques do not typically dominate migration 
(e.g., equation (l4TTi). 

But, recent work has suggested that nonisothermal be- 
havior could have an important effect on coorbital torques 
and the overall planet migration rate. The nonisothermal 
regime has been explored in simulations by Paardekooper 
& Mellema, 2006 who found that slowed and even outward 
migration due to coorbital torques may occur in certain 
regimes. Recent studies have suggested that a background 
radial entropy gradient could play a role in determining the 
corotational torque if the gas undergoes adiabatic changes 
in its interactions with the planet. This effect may provide 
an additional contribution to the coorbital torque beyond the 
vortensity gradient that appears in equation (|35] |. The slow- 
ing/reversing of migration appears to involve the conditions 
that the disk have a negative radial entropy gradient, suf- 
ficient viscosity to avoid coorbital torque saturation, and a 
thermal timescale that is long enough for the gas to behave 



adiabatically as it passes the planet. The nature of this ef- 
fect and the conditions required for it to operate is an active 
area of investigation (e.g., Baruteau & Masset, 2008). 

3.2 Other Migration Processes Involving a Gaseous 
Disk 

Runaway Coorbital Migration, Type III Migration 

The coorbital torque will be saturated (reduced to zero) 
unless some effect is introduced to break the symmetry of 
the streamlines, such as turbulent viscosity (see discussion 
in Section 2.2). However, planet migration itself introduces 
an asymmetry and could therefore act to prevent torque sat- 
uration. (The coorbital torque model presented in Section 
2.3 ignored the effects of planet migration on the horseshoe 
orbits in Fig.|2]) The coorbital torque for a migrating planet 
could then depend on the rate of migration. Under some 
conditions, the coorbital torque could in turn cause faster 
migration and in turn a stronger torque, resulting in an insta- 
bility and a fast mode of migration {Masset & Papaloizou, 
2003). The resulting migration is sometimes referred to as 
Type 111 migration. To see how this might operate in more 
detail, we consider the evolution of gas trapped in the coor- 
bital region {Artymowicz, 2004; Ogilvie & Lubow, 2006). 
For a sufficiently fast migrating planet, the topology of the 
streamlines changes with open streamlines flowing past the 
planet and closed streamlines containing trapped gas (see 
Fig. [T6] l. The leading side of the planet contains trapped 
gas acquired at larger radii, while the gas on the trailing 
side is ambient material at the local disk density. The den- 
sity contrast between material on the trailing and leading 
sides of the planet gives rise to a potentially strong torque. 
The migration timescales are of order the Type I migration 
timescales. This timescale can be very short for planets 
that are massive enough to partially open a gap and would 
otherwise not undergo Type I migration. The migration 
timescales have been found to be as short as 10-20 orbits 
in the case of a Saturn mass planet embedded in a cold and 
massive disk. A major question centers around the condi- 
tions required for this form of migration to be effective and 
therefore whether planets typically undergo such migration. 

Migration of Eccentric Orbit Planets 

The analysis of migration in Section 2 assumed that 
planets reside on circular orbits. This assumption is not 
unreasonable, since there are strong damping effects on 
eccentricity for a planet that does not open a gap in the 
disk (Artymowicz, 1993). Some eccentricity may be con- 
tinuously produced by turbulent fluctuations in the gas, as 
described above, or by interactions with other planets. In 
general, eccentricity damping is faster than Type I migra- 
tion. For planets that open a gap, it is possible that they re- 
side on eccentric orbits in the presence of the gaseous disk. 
In fact, one model for the observed orbital eccentricities of 
extra-solar planets attributes the excitation of eccentricities 
to disk-planet interactions {Goldreich & Sari, 2003; Ogilvie 
& Lubow, 2003). 

A planet on a sufficiently eccentric orbit embedded in a 
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circular disk can orbit more slowly at apoastron than the ex- 
terior gas with which it tidally interacts. Similarly, a planet 
can orbit more rapidly than the tidally interacting gas at 
periastron. These angular velocity differences can change 
the nature of the "friction" between the planet and the disk 
discussed in Section 2.2. For example, at apoastron the 
more slowly orbiting planet could gain angular momentum 
from the more rapidly rotating nearby gas that lies outside 
its orbit. Furthermore, since the planet spends more time 
at apoastron than periastron, the effects at apoastron could 
dominate over effects at periastron. It is then possible that 
outward migration could occur for eccentric orbit planets 
undergoing Type 1 migration, assuming such planets could 
maintain their eccentricities (Papaloizou, 2002). 

In the case of a planet that opens a gap, simulations sug- 
gest that slowed or even outward migration may occur as 
the planet gains eccentricity from disk-planet interactions 
(D' Angela et al, 2006). The situation is complicated by 
the fact that the gaseous disk can gain eccentricity from the 
planet by a tidal instability {Lubow, 1991). For the outward 
torque to be effective, there needs to be a sufficient differ- 
ence in the magnitude and/or orientation between the planet 
and disk eccentricities, so that the planet moves slower than 
nearby disk gas at apoastron. 

Multiplanet Migration 

Thus far, we have only considered single planet systems. 
Of the more than 200 planetary systems detected to date 
by Doppler techniques, over 10% are found in multi -planet 
systems {Butler et al, 2006). Systems have been found to 
have orbits that lie in mutual resonance, typically the 2:1 
resonance. The resonant configurations are likely to be the 
result of convergent migration, migration in which the sepa- 
ration of the orbital radii decreases in time. This process oc- 
curs as the outer planet migrates inward faster than the inner 
planet. Planets can become locked into resonant configura- 
tions and migrate together, maintaining the planetary orbital 
frequency ratio of the resonance (see also Chapter 10). The 
locking can be thought of as a result of trapping the plan- 
ets within a well of finite depth. Just which resonance the 
planets become locked into depends on their eccentricities 
and the relative rate of migration that would occur if they 
migrated independently. As planets that are initially well- 
separated come closer together, they lock into the first res- 
onance that provides a deep enough potential to trap them 
against the effects of their convergence. We discuss below 
some consequences of resonant migration. 

To maintain a circular orbit, a migrating planet must ex- 
perience energy, E, and angular momentum, J. changes 
that satisfy E ~ ftpj, where the angular speed of the planet 
rip varies as the planet migrates. As the planets migrate 
together in a resonant configuration, their mutual interac- 
tions cause deviations from this relation. As a result, their 
energies and angular momenta evolve in a way that is in- 
compatible with maintaining a circular orbit. Substantial 
orbital eccentricities can develop as a consequence of mi- 
gration (Yu & Tremaine, 2001; Lee & Peale, 2002). In addi- 
tion, migration can cause a large amplification of an initially 



small mutual inclination of the planetary orbit planes (Yu & 
Tremaine, 2001; Thommes & Lissauer, 2003; see Fig. fTTb . 

Planetary system GJ876 is a well-studied case in which 
the planets are in a 2:1 resonance. If the system's measured 
eccentricities are due to resonant migration, then according 
to theory (Lee & Peale, 2002), the system migrated inward 
by less than 10%. Such a small amount of locked migration 
seems unlikely. It is more reasonable to expect that some 
process limited further eccentricity growth. Disk-planet in- 
teractions could have limited the eccentricities and inclina- 
tions that could be developed by resonant migration. For a 
single planet interacting with a disk, theory and simulations 
suggest that eccentricity is generally damped for planets of 
low mass, too low to open a gap. A higher mass (gap open- 
ing) planet can undergo eccentricity growth due to its inter- 
action with the disk. But the level of eccentricity produced 
is limited and eccentricity damping occurs above that level 
{D' Angela et al, 2006). The disk-planet mutual inclination 
can also be suppressed or limited by the dissipation of disk 
warps (Lubaw & Ogilvie, 2001). Recent simulations sug- 
gest that disk-planet interactions due to gas interior to the 
orbit of the inner planet could have limited the eccentrici- 
ties of multi-planet systems to observed levels (Crida et al. 
2008). But it may be possible under certain circumstances, 
such as a depleted inner disk, that large eccentricities and 
orbital inclinations could develop due to such processes. 

3.3 Other Forms of Migration 

Migration in a Planetesimal Disk 

After the gaseous disk is cleared from the vicinity of the 
star, after about lO'^y, there remains a disk of solid material 
in the form of low mass planetesimals. This disk is of much 
lower mass than the original gaseous disk. But the disk is 
believed to have caused some migration in the early solar 
system with important consequences (Hahn & Malhotra, 
1999; Tsiganis et al, 2005). 

There is strong evidence that Neptune migrated outward. 
This evidence comes from observations of Kuiper belt ob- 
jects that are resonantly trapped exterior, but not interior, 
to Neptune's orbit. The detailed dynamics of a planetesi- 
mal disk are somewhat different from the case of a gaseous 
disk, as considered in Section 2. The planetesimals behave 
as a nearly collisionless system of particles. Jupiter is much 
more massive than the other planets and can easily absorb 
angular momentum changes in Neptune. As Neptune scat- 
ters planetesimals inward and outward, it undergoes angular 
momentum changes. It is the presence of Jupiter that breaks 
the symmetry in Neptune's angular momentum changes. 
Once an inward scattered planetesimal reaches the orbit of 
Jupiter, it gets flung out with considerable energy and does 
not interact again with Neptune. As a result of the loss of in- 
ward scattered particles, Neptune gains angular momentum 
and migrates outward, while Jupiter loses angular momen- 
tum and migrates slightly inward. 

A somewhat analogous process occurs in gaseous cir- 
cumbinary disks, disks that orbit around binary star systems 
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(Pringle, 1991). The circumbinary disk gains angular mo- 
mentum at the expense of the binary. The binary orbit con- 
tracts as the disk outwardly expands. A gap-opening planet 
embedded in a circumbinary disk (or under some condi- 
tions, a disk that surrounds a star and massive inner planet) 
would undergo a form of Type II migration that could carry 
the planet outward {Martin et al, 2007). In the solar sys- 
tem case, the Sun- Jupiter system plays the role of the bi- 
nary. Viscous torques are the agent for transferring the an- 
gular momentum from the binary outward in the gaseous 
circumbinary disk, while particle torques play the some- 
what analogous role in the planetesimal disk. 

In a planetesimal disk, another process can operate to 
cause migration. This process is similar to the runaway 
coorbital migration (Type III migration) described above, 
but applied to a colUsionless system of particles (Ida et 
al. 2000). Interactions between the planetesimals and the 
planet in the planet's coorbital zone can give rise to a mi- 
gration instability. 

Kozai Migration 

A planet that orbits a star in a binary star system can 
periodically undergo a temporary large increase in its or- 
bital eccentricity through the a process known as the Kozai 
effect (see also Chapter 10). Similar Kozai cycles can oc- 
cur in multi-giant planet systems. The basic idea behind 
Kozai migration is that the increased eccentricity brings the 
planet closer to the star where it loses orbital energy through 
tidal dissipation. In the process, the planet's semi-major 
axis is reduced and inward migration occurs (Wu & Mur- 
ray, 2003). We describe this in more detail below. 

Consider a planet in a low eccentricity orbit that is 
well interior to the binary orbit and is initially highly in- 
chned with respect to it. The orbital plane of the planet 
can be shown to undergo tilt oscillations on timescale of 
^ P^/Pp, where Pb is the binary orbital period, Pp is the 
planet's orbital period, and by assumption P^ Pp. Un- 
der such conditions, it can be shown that the component 
of the planet's angular momentum perpendicular to the bi- 
nary orbit plane (the z-component) is approximately con- 
served, Jz = Mp^GMsap{l — Cp) cos/, where ap and 

Bp are respectively the semi-major axis and eccentricity of 
the planet's orbit and / is the inclination of the orbit with 
respect to the plane of the binary. 

The conservation of is easily seen in the case that 
the binary orbit is circular and the companion star is of 
low mass compared to the mass of the star about which the 
planet orbits. On such long timescales Pb> the compan- 
ion star can be considered to be a continuous ring that pro- 
vides a static potential. In that case, the azimuthal symme- 
try of the binary potential guarantees that is conserved. 
By assumption, we have cos / <C 1 and gp ^ 1 in the ini- 
tial state of the system. As the planet's orbital plane evolves 
and passes into aUgnment with the binary orbital plane, 
cos/ ^ 1, conservation of Jz requires Cp ^ 1. In other 
words, Jz is initially small because of the high inclination 
of the orbit. When the inclination drops, the orbit must be- 



come more eccentric (radial), in order to maintain the same 
small Jz value. The process then periodically trades high 
inclination for high eccentricity. 

During the times of increased eccentricity, the planet 
may undergo a close encounter with the central star at peri- 
astron distance ap(l — ep). During the encounter, the tidal 
dissipation involving the star and planet results in an energy 
loss in the orbit of the planet and therefore a decrease in Op. 
This process then results in inward planet migration. The 
energy loss may occur over several oscillations of the orbit 
plane. 

Another requirement for the Kozai process to operate is 
that the system must be fairly clean of other bodies. The 
presence of another object could induce a precession that 
washes out the Kozai effect. Given the special requirements 
needed for this process to operate, it is not considered to 
be the most common form of migration. However, there is 
good evidence that it does operate in some systems (Takeda 
& Rasio, 2005). The Kozai effect can also occur in two- 
planet systems, where the outer planet plays the role of the 
binary companion. The process can be robust due to the 
proximity of outer planet (Nagasawa et al. 2008). 

3.4 Techniques Used in Numerical Simulations 

Numerical simulations provide an important tool for an- 
alyzing planet migration. They can provide important in- 
sights in cases where nonlinear and time-dependent effects 
are difficult to analyze by analytic methods. Simulations, 
as well as analytical models, depend on a model of appro- 
priate physical processes, such as heating and cooling, the 
treatment of the disk turbulence (either by a viscosity or de- 
tailed modehng of the disk instability that causes the turbu- 
lence), and the model of the disk structure. Below we dis- 
cuss some issues related to the use of simulations for sim- 
ple (isothermal and alpha-disk) multi-dimensional models 
of disk-planet interactions. 

Some powerful grid-based hydrodynamics codes (such 
as the ZEUS code Stone & Norman, 1992) have been 
adapted to the study of disk-planet interactions. In addi- 
tion, particle codes based on the Smoothed Particle Hydro- 
dynamics (SPH) (Monaghan, 1992) have sometimes been 
employed. A systematic comparison between many of the 
codes has been carried out by de Val-Borro et al, 2006. We 
discuss a few basic points. For planets that open a gap, grid- 
based codes offer an advantage over particle-based codes. 
The reason is that the resolution of grid based codes is deter- 
mined by the grid spacing, while the resolution of particle- 
based codes is determined by the particle density. If a planet 
opens an imperfect gap, the particle density and resolution 
near the planet is low. Higher resolution occurs where the 
particle density is higher. On the other hand, SPH well rep- 
resents regions of a cold disk that are away from a planet. 
The gas follows the trajectories of SPH particles that are 
slightly modified by gas pressure and tidal forces. SPH 
provides high resolution in dense regions, such as near the 
cores of low mass planets embedded in disks (Ayliffe & 
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Bate, 2007). 

There have been variable resolution techniques devel- 
oped for grid-based codes in which the highest resolution 
is provided in regions near the planet where it is needed, 
such as nested grid methods (D' Angela et al, 2002). Such 
techniques need to provide a means of joining the regions of 
high and low resolution without introducing artifacts (such 
as wave reflections) or lowering the overall accuracy of the 
scheme. 

Grid-based codes that simulate disk-planet interactions 
typically employ numerical devices to improve conver- 
gence. For example, the gravitational potential of the planet 
is often replaced by one that does not diverge near the 
planet. The potential is limited by introducing a smooth- 
ing length, a distance within which the potential does not 
increase near the planet. Another limitation is that the sim- 
ulated domain of the disk is typically limited to a region 
much smaller than the full extent of the disk. Techniques 
have been developed to ensure that reflections from the 
boundaries do not occur, e.g., by introducing enhanced 
wave dissipation near the boundary or approximate out- 
going wave boundary conditions. 

The time-steps of codes are limited by the Courant con- 
dition. Short time-steps often result from the region near 
the inner boundary of the computational domain (smaller 
radii) where the disk rotation is fastest. As a result, it is dif- 
ficult to extend the disk very close to the central star The 
FARGO scheme (Masset, 2000) is very useful in overcom- 
ing this limitation. However, the method is difficult to apply 
to a variable grid spaced code. 

Convergence is a major issue with these simulations. 
Ideally, one should demonstrate that the results of the sim- 
ulations are sufficiently insensitive to the locations of the 
boundaries, the size of the smoothing length, the size of 
the time steps, and the grid resolution. Even the direction 
of migration can be affected by the size of the potential 
smoothing length in certain cases. In a 2D simulation, a 
finite smoothing length ^ H provides a means of simulat- 
ing the reduced effects of planet gravity on a disk of finite 
thickness. But in the 2D case, the limit of zero smoothing 
length is unphysical. In practice, testing for convergence is 
computationally expensive, but can be done for a subset of 
the models of interest, perhaps over a limited time range. 
Demonstrating convergence is important for providing reli- 
able results. 



4. FUTURE PROSPECTS 
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Fig. 1. — Migration timescales versus planet mass for a planet 
embedded in a 3D disk of mass ~ 0.02Ma with surface density 
E oc r~^^^ and H/r = 0.05. The planets are on fixed circular 
orbits and have fixed masses and orbit a solar mass star. The dots 
with error bars denote results of 3D numerical simulations with 
the same disk parameters {Bate et al, 2003). The dashed line plots 
equation i4H based on linear theory (Tanaka et al, 2002). Above 
about O.lAfj, the planet opens a gap in the disk. Type I theory 
becomes invalid, and Type II migration occurs. 
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The theory of planet migration is intertwined with the 
theory of planet formation (Chapter 13). The timescale for 
a planet to grow within a disk is a key element in under- 
standing whether planet migration is a major obstacle to 
planet formation. We have pointed out that the formation 
timescales for gas giant planets in the core accretion model 
do present a problem for the simplest planet migration theo- 
ries. However, alternative migration models, some of which 



Fig. 2. — Coorbital streamlines near a Saturn mass planet that 
orbits a solar mass star, Mp = 3 x 10~''Ms. The origin is centered 
on the planet and corotates with it. Azimuthal angle 9 about the 
star increases to the left. The outermost streamlines are at the edge 
of the coorbital region. Solid streamlines are the horseshoe orbits. 
Locations i=p and o=p label the positions near the encounter with 
the planet. 
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Fig. 3. — Path of a particle that passes by a planet of mass Afp = 
10"'^Ms (0.3 Eai-th masses for a planet that orbits a solar mass 
star). The coordinates are in units of the orbital radius of the planet 
a. The planet lies at the origin, while the star lies at ( — 1,0). The 
dashed line follows the path that is undisturbed by the planet with 
X = 0.025a, while the solid line follows the path resulting from 
the interaction with the planet. In the frame of the planet, the 
particle moves in the negative y direction. 




Fig. 4. — The solid line plots the particle distance from the star r 
in units of a as a function of time in units of planet orbit periods 
for the particle that follows the perturbed path in Fig.|3] The par- 
ticle passes the planet at time t = Q. Immediately after passage 
by the planet, the particle is deflected toward smaller radii, toward 
the planet, and acquires an eccentricity, as indicated by the radial 
oscillations. The dashed line plots the mean radius of these os- 
cillations. Since the mean radius of the oscillations is larger than 
the initial orbital radius (compare dashed line with solid line at 
t < — 1), the particle gained energy and angular momentum, as a 
consequence of its interaction with the planet. 



are described in Section 3, may be appropriate. Future 
prospects for resolving this issue depend on advances in the 
theory of planet formation. 

Prospects for making progress in the theory of planet mi- 
gration rely on including and more accurately representing 
physical processes in the models, improving computer sim- 
ulations, and obtaining a better understanding of the struc- 
ture of planet-forming gaseous disks. 

Planet migration in disks is a consequence of the action 
of Lindblad and corotational resonances, as described in 
Section 2. The theory of Lindblad resonances is better un- 
derstood. Their linear and nonlinear properties have been 
analyzed in more detail. The corotational resonances are 
somewhat more delicate and less well understood. Progress 
on the theory of planet migration will likely involve further 
investigations of the role of corotational resonances (e.g., 
Paardekooper & Mellema, 2006). 

Improvements to computer capabilities and codes should 
allow simulations to be caiTied out with higher resolution 
over longer timescales. Progress will be made by also in- 
cluding more physical effects. For example, most multi- 
dimensional simulations have made only the simplest as- 
sumptions about the thermal properties of the disk. The use 
of a turbulent viscosity is a serious approximation to the ef- 
fects of turbulence. Directly simulating the instability that 
produces the turbulence along with planet migration is a 
major computational challenge (e.g.. Nelson, 2005). 

Some calculations have suggested that higher mass ec- 
centric orbit planets could undergo slowed or outward mi- 
gration due to their interactions with a gaseous disk. But 
we do not know whether such planets have acquired their 
eccentricities at this early stage. Observations of the ec- 
centricities young planets would be quite valuable in un- 
derstanding this issue. It would be useful to know whether 
planetary eccentricities are determined at early times, sug- 
gesting that planet eccentricities are present while the planet 
is embedded in its gaseous disk. 

A major uncertainty in the theory of planet migration is 
the physical state of the disk. Low mass planet migration 
behaves very differently depending on the level of disk tur- 
bulence and the structural properties of the disk (Section 3). 
For example, in weakly turbulent disks, feedback effects 
may halt Type I migration. While in highly turbulent disks, 
fluctuating torques may dominate over Type I torques. The 
presence of rapid radial density variations or weakened (or 
inverted) radial temperature variations can substantially al- 
ter migration, since it depends on the competition between 
torques involving material just inside and outside the orbit 
of the planet. The better determination of disk properties 
will likely rely on some combination of improved theory 
and observations. It is unlikely that theory alone will be 
able to make much progress along these lines because, for 
example, the density structure depends on the disk turbu- 
lence that is difficult to accurately predict from first princi- 
ples. Once such difficulty is that magnetically driven turbu- 
lence is strongly affected by the abundance and size distri- 
bution of dust grains that control the level of disk ionization 
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(e.g., Salmeron & Wardle, 2008). The observational deter- 
mination of disk properties (disk density and temperature 
as a function of distance from the star) in the inner parts 
of protostellar disks, where planet formation is expected to 
occur, is important for such purposes. New telescopes such 
as ALMA and JWST may be quite valuable in making such 
determinations. 
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Fig. 5. — Torque on the particle of mass M that is normalized 
by ICPMQ.I a? as a function of time in units of the planet's orbit 
period along the unperturbed (dashed lines) and perturbed paths 
(solid lines). The planet mass Mp = 10~^Afs (0.3 Earth masses 
for a planet that orbits a solar mass star). Top panel is for a particle 
having x — 0.02a, the case in Fig. |3] The bottom panel is for a 
particle with x = 0.03a. 




Fig. 6. — Numerical test of equation ((S} based on orbit integra- 
tions. Top panel: Log-log plot of 10*A J/ (Ala^flp) as a function 
of preimpact distance from the planet's orbit x in units of a. The 
lower set of points is for a planet of mass Mp — 10~^Afs (0.3 
Earth masses for a planet that orbits a solar mass star). The upper 
set is for Alp = 10^'^A/s (a Jupiter mass for a planet that orbits 
a solar mass star). The solid lines are for A J oc x^'^ that pass 
through the respective right-most points. Bottom panel: Log-log 
plot of 10'' A J/(Ma^r2p) as a function of planet mass in Earth 
masses for a planet that orbits a solar mass star. The points are 
the results of numerical simulations for a fixed value of a; = 0.2o. 
The solid line is for AJ oc (Mp/Ms)^ that passes through the 
left-most point. 
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Fig. 7. — Total angular momentum for a set of 60 particles, each 
of mass M, on horsesiioe orbits in units of Ma?Q,p as a function 
of time in units of the planet orbit period in a star-planet system 
with Mp = 3 X 10~®Ms (0.3 Earth masses for a planet that orbits 
a solar mass star). The particles start at t = distributed between 
r = a + 7?h/60 to r = a + iln (with _Rh ^ 7 x lO^^a) along 
the star-planet axis 180 degrees from the planet. The changes in 
angular momentum cause a torque to be exerted on the planet. The 
torque (the time derivative of J) oscillates and declines in time as 
angular momentum becomes more constant in time because the 
particles undergo phase mixing on the libration timescale ~ 150 
planet periods. 
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Fig. 9. — Schematic of acoustic (pressure) wave propagation 
in a gas disk involving an inner Lindblad resonance (ILR), outer 
Lindblad resonance (OLR), and corotation resonance (CR). Spi- 
ral waves launched at Lindblad resonances propagate away from 
the orbit of the planet. The region near the CR, between ILR and 
OLR, is evanescent (nonpropagating). 
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Fig. 8. — Radius of a Lindblad resonance in units of a as a 
function of azimuthal wave number m. Points below (above) the 
dashed line are for inner (outer) Lindblad resonances. The dashed 
line is the location of the planet. The two sets of smaller dots 
are for a Keplerian disk with radii given by equation l |38l ). The 
resonances get closer to the planet with increasing m. The two 
sets of larger dots account for azimuthal pressure effects in a disk 
with thickness to radius ratio H /r = 0.1. With these pressure 
effects included, the resonances maintain a fixed separation from 
the planet at high m. 
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Fig. 10. — Migration of a planet orbiting a solar mass star and 
undergoing growth via gas accretion. The disk parameters are sim- 
ilar to those in Fig[T] Top: The vertical axis is the orbital radius 
in units of the initial orbital radius 5.2AU. The horizontal axis is 
time in units of the initial orbital period 12yr. The solid curve is 
the result of 3D hydrodynamical simulations. The lower and upper 
short dashed curves are based on equations HOt and l l41t respec- 
tively, applied to a planet of variable mass. The long dashed curve 
corresponds to migration on the disk viscous timescale. Bottom: 
Average disk density near the planet relative to the initial value 
as a function of time. The density is averaged over a band of 
radial width 2H centered on the orbit of the planet and is nor- 
malized by its initial value. The first solid circle marks the time 
when Mp = 16.7M^, subsequent circles occur at integer mul- 
tiples of 33M0. The planet initially follows the predictions of 
Type I migration theory in the top panel while there is no substan- 
tial gap in the disk (no drop in the curve on the bottom figure). 
After gap opening, the planet follows Type II migration. Obtained 
from D' Angela & Lubow, 2008. 




Fig. 11. — Scaled torque per unit disk mass on the planet as a 
function of radial distance from the planet's orbital radius based 
on 3D simulations. The horizontal scale is in units of a and the 
vertical scale is in units of GMs{M-p/Ms)'^ /a. The solid and 
long-dashed curves are for 1 (1 Earth mass) and 10 mass 
planets, respectively, that orbit a solar mass star. The disk param- 
eters are H/r = 0.05 and a = 0.004 for both the cases. Ac- 
cording to linear theory, these two curves should overlap. Torque 
distributions are averaged over one orbital period. Figure based on 
D' Angela & Liibaw, 2008. 
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Fig. 12. — Azimuthally averaged disk surface density normalized 
by the unperturbed value at radius a as a function of r in units of 
a = 5.2AU. The central star has a solar mass. The disk is sim- 
ulated in 3D with parameters H/r — 0.05 and a = 0.004. The 
density profiles are for planets with masses of 1 (long-dashed), 
0.3 (dot-dashed), 0.1 (dotted), 0.03 (short-dashed), and 0.01 (thin 
solid) A/j. Only planets with masses Mp > 0.1 A/j produce sig- 
nificant perturbations. The thick solid line is based on a 2D simu- 
lation of a lA'/j planet by Lubow et al (1999). Obtained from Bate 
et al, 2003. 
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Fig. 13. — Schematic of the axisymmetric density perturbation 
in arbitrary units as a function of radial distance from the planet 
(shown as dot at the origin) in units of disk thickness H. The 
solid lines with arrows indicate the velocity perturbation on the 
gas caused by tidal torques on the planet. These velocities are 
directed away from the orbit of the planet. The dashed lines with 
arrows indicate the gas velocity, in the frame of the planet, due to 
the assumed inward migration of the planet. Interior to the orbit of 
the planet, the velocities due to tides and migration oppose each 
other and lead to a pile up of material. Outside the orbit of the 
planet, the velocities reinforce each other and lead to a density 
drop. 
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Fig. 15. — Orbital radius as a function of time in orbits for SM^ 
planets embedded in a gaseous disk based on simulations. The 
dotted lines plot the migration of planets in a disk without turbu- 
lent fluctuations. The solid lines plot the migration of planets in a 
disk with turbulent fluctuations due to the MHD turbulence. The 
random motions are due to the fluctuating torques. Obtained from 
Nelson, 2005. 



0.75 



a =8 X 10" 




200 400 600 800 1000 1200 

Orbits 



Fig. 14. — Influence of disk viscosity parameter a on the migra- 
tion of a planet with mass lOAf0 in a disk with Hjr = 0.035 and 
mass ~ O.IMq based on 2D simulations. The vertical axis is the 
orbital radius in units of the initial orbital radius. The horizontal 
axis is the time in units of the initial planet orbital period. For 
a = 8 X 10~^ the migration follows the Type 1 rate. At the lower 
values, the migration halts due to a feedback effect (see Fig. [T3] l. 
Figure based on Li et al, 2008. 
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Fig. 16. — Coorbital streamlines near an inwardly migrating Sat- 
urn mass planet Mp = 3 x 10~*A/s located at the origin in the 
comoving frame of the planet that orbits a solar mass star. The 
inward migration rate is 0.002f2pa. Angle 9 increases to the left. 
The heavy streamlines are open and pass by the planet. The light 
streamlines are closed and contain trapped material on the leading 
side of the planet (streamlines based on Ogilvie & Lubow, 2006). 
The trapped material, acquired at earlier time, is retained from re- 
gions further from the star. The open streamlines carry ambient 
disk material. In contrast to the nonmigrating case in Fig. (2] the 
asymmetry and the density differences between the trapped (re- 
tained) and open (ambient) gas gives rise to a coorbital torque. 
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Fig. 17. — Evolution of orbital eccentricity (top) and inclination 
(bottom) of a two-planet system in which an inward migration of 
10~^AUyr^^ is imposed on the outer one. The horizontal axis is 
in units of years. Both planets have a mass equal to Jupiter's mass. 
The orbital eccentricities of the planets grow as a result of resonant 
migration in the 2: 1 resonance. The eccentricity of the inner planet 
grows faster than the outer one. A small initial mutual inclination 
grows to large values as the eccentricities increase. From Tommes 
& Lissauer, 2003. 
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